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ABSTRACT 

The European Gaia astrometry mission is due for launch in 2011. Gaia will rely on the 
proven principles of ESA’s Hipparcos mission to create an all-sky survey of about one billion 
stars throughout our Galaxy and beyond, by observing all objects down to 20th magnitude. 
Through its massive measurement of stellar distances, motions and multi-colour photometry it 
will provide fundamental data necessary for unravelling the structure, formation and evolution 
of the Galaxy. This paper presents the design and performance of the broad- and medium-band 
set of photometric filters adopted as the baseline for Gaia. The nineteen selected passbands 
(extending from the ultraviolet to the far-red), the criteria, and the methodology on which 
this choice has been based are discussed in detail. We analyse the photometric capabilities 
for characterizing the luminosity, temperature, gravity and chemical composition of stars. We 
also discuss the automatic determination of these physical parameters for the large number of 
observations involved, for objects located throughout the entire Hertzsprung-Russell diagram. 
Finally, the capability of the photometric system to deal with the main Gaia science case is 
outlined. 

Key words: Galaxy: structure, formation - stars: fundamental parameters - techniques: pho¬ 
tometric - instrumentation: photometers - space vehicles: instruments 


2 Jordi et al 


1 INTRODUCTION 

Gaia has been approved as a cornerstone mission in the ESA sci¬ 
entific programme. The main goal is to provide data to study the 
formation and subsequent dynamical, che mical and star f ormation 
evolu tion of the Milky Way galaxy j Perry man et al 1200 it [m i ” n a rd 
l 206 l . Gaia will achieve this by providing an all-sky astrometric 
and photometric survey complete to 20 mag in unfiltered light. 
During the mission, on-board object detection will be employed 
and more than 1 billion stars will be observed (as well as non- 
stellar objects to similar completeness limits). The full mission 
(5-years) mean-sky parallax accuracies are expected to be around 
7 microarcsec at V = 10, 12-25 microarcsec at V = 15 and 100- 
300 microarcsec at V = 20 (depending on spectral type). Multi¬ 
epoch, multi-colour photometry covering the optical wavelength 
range will reach the same completeness limit. Radial velocities will 
be obtained for 100-150 million stars brighter than V ~ 17—18 mag 
with accuracies of around 1-15 km s -1 , depending on the apparent 
magnitude and s pectral type of the stars and the sky density (for 
details see lKatz et all2004l : IWilkinson et all2005l) . 

The photometric measurements provide the basic diagnos¬ 
tics for classifying all objects as stars, quasars, solar-system ob¬ 
jects, or otherwise, and for parametrizing them according to their 
nature. Stellar classification and parametrization across the entire 
Hertzsprung-Russell diagram is required as well as the identifica¬ 
tion of peculiar objects. This demands observation in a wide wave¬ 
length range, extending from the ultraviolet to the far-red. The pho¬ 
tometric data must determine: (i) effective temperatures and red¬ 
dening at least for O-B-A stars (needed both as tracers of Galactic 
spiral arms and as reddening probes), (ii) at least effective temper¬ 
atures and abundances for F-G-K-M giants and dwarfs, (iii) lumi¬ 
nosities (gravities) for stars having large relative parallax errors, 
(iv) indications of unresolved multiplicity and peculiarity, and (v) 
a map of the interstellar extinction in the Galaxy. All of this has to 
be done with an accuracy sufficient for stellar age determination in 
order to allow for a quantitative description of the chemical and dy¬ 
namical evolution of the Galaxy over all galactocentric distances. 
Separate determination of Fe- and a-element abundances is essen¬ 
tial for mapping Galactic chemical evolution and understanding the 
formation of the Galaxy. 

Photometry is also crucial to identify and characterize the set 
of ~ 500 000 quasars that the mission will detect. Apart from be¬ 
ing astrophysically interesting in their own right, quasars are key- 
objects for defining the fixed, non-rotating Gaia Celestial Referen¬ 
ce Frame, the optical equivalent of the International Celestial Refer¬ 
ence Frame fMignartfcoQsth On the other hand, Gaia will identify 
about 900 quasars with multiple images produced by macrolens- 
ing. Since this number is sensitive to cosmological parameters, the 
Gaia observations will be able to constrain the latter. 

Due to diffraction and the optical aberrations of the instru¬ 
ment, the position of the centre of the stellar images is wavelength 
dependent. To achieve the microarcsec accuracy level, astrometry 
has to be corrected for this chromatic aberration through the knowl¬ 
edge of the spectral energy distribution of the observed objects. 
Photometry is indispensable for this. If uncorrected, chromatic er¬ 
rors could reach several milliarcsec, cf. Section l5.1.2l 

As explained in the following sections, the photometric sys¬ 
tems proposed during the long development of Gaia have been 
improved along with the increasing collecting area of the tele¬ 
scopes, with better insight in the astrophysical requirements, and 
with the development of mathematical tools to compare the vari¬ 
ous proposed systems. The use of CCDs in a scanning astrome¬ 


try satellite was first proposed in 1992 fae)d ll993) as the ROH¬ 
MER project. The proposal included five broad passbands, UB- 
VRI, which would obtain much better precision than the B, V of 
Hipparcos-Tycho although with a similar collecting aperture. The 
Gaia collecting area has increased by up to 10 times with respect to 
ROEMER and consequently the initial photometric system has been 
upgraded severa l times . Eight m edium-width passbands were pro¬ 
posed by Straizv s~& Hodi 19951) . and sp ectrophotom etry instead of 
filter photometry was also considered in H^g| 1 1 9981) . A system of 
four broad and elev en medium-width passbands was propos ed by 
[Grenon et ah ( 1999 ) and adopted in the ‘Gaia Study Report’ l lESAl 
|200C|). The subsequent developments and updates have yielded the 
present baseline with five broad- and fourteen medium-width pass- 
bands. 

This paper deals with the definition of the photometric sys¬ 
tem, the relationship of its passbands with the stellar astrophysical 
diagnostics and the evaluation of its performance in terms of the 
astrophysical parametrization of single stars. This is based on the 
specific design implementation of the payload commonly referred 
to as Gaia-2 (cf. Section|3 applicable at the end of the technology 
assessment phase as of mid-2005. The resulting astrometric and 
photometric requirements form the basis of the industrial specifica¬ 
tions for the satellite implementation phase, with the consequence 
that the detailed design, due for finalisation early in 2007, may dif¬ 
fer in detail from the present description. Nevertheless, the prin¬ 
ciples and objectives as well as the methods and assessment tools 
described in this paper will remain applicable. 

The paper is organized as follows: Section[2]describes the mis¬ 
sion observation strategy, telescopes and focal planes. Section [3] 
deals with the measurement of the unfiltered light in the fields of 
view. The principles of designing the multi-colour photometric sys¬ 
tem are outlined in Section [4] and the purpose of the broad- and 
medium-passbands is discussed in detail in Section [5] Synthetic 
photometry and corresponding error estimates are given in Sec- 
tion[6] The performance of the photometric system with respect to 
astrophysical parameter determination is quantified in Section^ In 
Section|8|the potential of the photometric system for Galactic struc¬ 
ture and evolution studies as well as the performances for QSOs 
are outlined. Finally, Section [9] and Appendix [A] present the con¬ 
clusions and describe the ‘Figure of Merit’ mathematical tool for 
the comparison of filter systems, respectively. 

Throughout the rest of the paper we will use the following 
abbreviations: AP: astrophysical parameter; BBP: broad-band pho¬ 
tometer; C1B: set of passbands implemented in the Astro instru¬ 
ment; C1M: set of passbands implemented in the Spectra instru¬ 
ment; CCD: charge coupled device; FoM: figure of merit; MBP: 
medium-band photometer; MDM: minimum distance method; NN: 
neural network; PS: photometric system; QE: CCD quantum ef¬ 
ficiency; QSO: quasar; SED: spectral energy distribution; SNR: 
signal-to-noise ratio; ST: scientific target. 

2 INSTRUMENT DESCRIPTION 
2.1 Observation strategy 

Gaia is a survey mission. Operating from a Lissajous orbit around 
the second Lagrange point of the Sun-Earth/Moon system (L2), the 
satellite will continuously observe the sky. During its 5-year mis¬ 
sion, Gaia will rotate, at a fixed speed of 60 arcsec s -1 , around a 
slowly precessing spin axis. As a result of this spin motion, objects 
traverse the focal planes, which have viewing directions oriented 
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Figure 1. Schematic layout for the Astro field (left) and the Spectro/MBP field (right) in Gaia- 2. The designations BBP1-BBP4 and MBP01-MBP20 refer to 
the physical locations of the CCDs in the along-scan direction within the field, independent of their functionality (e.g., for detection or photometry) and the 
assignment of filters between them. The Spectro sky mapper, operating without filters, is physically located at MBP01, 02, 06, and 07. ASM and AF are the 
Astro sky mapper and the Astro field, respectively. 


Table 1. Gaia -2 characteristics of the Astro and Spectro instruments. 



Instrument 

Photometer 

Astro 

Broad-band 

Spectro 

Medium-band (blue) 

Medium-band (red) 

F 

Focal length [m] 

46.7 

2.3 

2.3 

A 

Telescope pupil area [m 2 ] 

1.4x0.5 = 0.7 

0.5x0.5 = 0.25 

0.5x0.5 = 0.25 

^3-sup 

Number of superimposed fields 

2 

1 

1 

To (A) 

Telescope transmittance^^ 

Ti s w 

0.8 x !*(A) 

TLW 

T 

Integration time [s] per single CCD 

3.31 

12.0 

12.0 

Q( A) 

CCD quantum efficiency 

CCD-green 

CCD-blue 

CCD-red 

r 

Total detection noise [e — sample -1 ] 

6.6 

7.3 

7.3 


Available number of filter strips 

4 (10 rows) 

10 (2 rows) 

6 (2 rows) 


Field of view across scan [deg] 

0.7236 

1.4685 

1.4685 

n obs 

Mean number of observations per strip (5 years)( 2 ) 

83 

84 

84 


Sample size on sky [arcsec x arcsec] 

0.0442x1.5912 

0.897x1.346 

0.897x1.346 


Double-star resolution [arcsec] 

0.05-0.1 

0.5-1 

0.5-1 


No-filter magnitude 

G 

- 

GS 


W 7)7 (A) denotes n reflections of material X. The factor 0.8 in MBP-blue comes from the dichroic beam splitter. 

The number of observations with Astro ranges from '-40 to 250 and in Spectro from '-50 to 220. 

(: 0 Sample size on the sky: A sample is a binning of CCD pixels in an area along-scan X across-scan. 


at right angles to the spin axis. The along-scan direction is defined 
to be the direction in which the images move in the focal plane, 
and the across-scan direction as the perpendicular direction. Thus, 
an object crosses the field of view in the along-scan direction at a 
roughly constant across-scan coordinate. 

Gaia’s observations are made with high-quality, large-format 
charge coupled devices (CCDs). These detectors are operated in 
time delay and integration mode, with charge images being trans¬ 
ported (clocked) in synchrony with optical images moving across 
the field due to the rotation of the satellite. Each CCD may, in prin¬ 
ciple, be equipped with a filter (an interference, standard coloured 
glass or multi-layer filter) defining, together with the telescope 
transmission and the CCD quantum efficiency (QE), a certain pho¬ 
tometric passband. By design, the arrangement of CCDs and filters 
in the focal planes ensures redundancy in case of failure. 

Stars brighter than ~12 mag pose a special challenge. Pixel 
saturation may be avoided for such objects by the activation of gates 
on the CCD, effectively reducing the CCD integration time. 


2.2 Gaia instruments 

Gaia will perform measurements by means of two physically dis¬ 
tinct instruments, with different viewing directions: the Astro in¬ 
strument, designed for astrometric and broad-band photometric ob¬ 
servations, and the Spectro instrument, used for medium-band pho¬ 
tometry and radial-velocity measurements. Astro and Spectro dif¬ 
fer in spatial resolution, available integration time, the number and 
type of passbands that can be used, and in telescope transmission 
and detector characteristics (Tableland Figure^. 

The astrometric focal plane incorporates three functions: (i) 
the sky mapper; (ii) the main astrometric field; and (iii) the broad¬ 
band photometer (BBP; Section l5~Tl . BBP is mainly aimed at sam¬ 
pling the spectral energy distribution of objects over a wide wave¬ 
length range to allow on-ground correction of image centroids mea¬ 
sured in the main astrometric field for systematic chromatic shifts 
caused by aberrations. In addition, BBP measurements contribute 
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to the astrophysical characterisation of objects, especially in dense 
stellar fields. 

Through implementation of a dichroic beam splitter, the Spec¬ 
tra instrument serves two distinct focal planes: one for the radial- 
velocity spectrometer (RVS), and one for the medium-band pho¬ 
tometer (MBP; Section l5*2l . The Spectro/MBP focal plane incor¬ 
porates two functions: (i) a sky mapper; and (ii) the MBP instru¬ 
ment. The main goal of MBP is to determine the astrophysical 
parameters of objects, which, in combination with the astrometric 
measurements, will enable astronomers to fulfil Gaia s main sci¬ 
ence objective. RVS observes a spectral region around 860 nm at 
a nominal resolution of 11 500. Its main aim is to determine ra¬ 
dial velocities for bright stars, down to V — 17-18 mag. For the 
brightest ones, astrophysical parametrization is also foreseen. 

The angular resolution of the instruments allows photometry 
to be obtained at stellar densities up to 750 000 stars/sq deg in the 
BBP and up to 200 000-400 000 stars/sq deg in the MBP (cf. Sec- 
tion 16.3> . The resolution limit for double stars is about 0.05-0.1 
arcsec in the Astro field and 0.5-1 arcsec in the MBP. 


2.3 Sky mappers 

Both sky mappers work in unfiltered light and allow autonomous 
object detection and confirmation, including rejection of prompt 
particle events (cosmic rays, solar protons, etc.). Detection and 
confirmation probabilities are a function of magnitude and object 
density in the field of view. They are effectively unity up to the 
survey limit (20 mag) dropping quickly to zero for fainter objects 
dArenou et a 1I2005H) . On-board object detection has several advan¬ 
tages: (i) sampling of detected obj ects can be limited to ‘windows', 
i.e., areas centred on the object iH0 a200511 : this allows to flush use¬ 
less pixels containing empty sky, which is of benefit to both CCD 
read-out noise and the data volume to be transmitted to ground; (ii) 
it allows the unbiased detection of all objects, assuring that the re¬ 
sulting catalogue will be complete to the survey limit; and (iii) un¬ 
predictable ‘peculiar objects', such as supernovae or solar-system 
objects, will be observed, if brighter than the detection limit. More¬ 
over, on-board detection is mandatory as no input catalogue exists 
which is complete to the survey limit at the Gaia spatial resolution. 


2.4 Astro/BBP 

Gaia’s astrometric observations are made without a filter in order 
to minimize photon noise. The mirror coatings and CCD QE effec¬ 
tively define a broad (white-light) passband, called G (Section^- 
The Astro focal plane receives the light from the superposition of 
two viewing directions on the sky, separated by the so-called ‘basic 
angle’, of order 100°. This superposition of the two fields of view 
will lead to a number of complications in the data processing, but 
throughout this paper we only consider the first order effect of the 
doubling of the diffuse sky-background. 

The BBP consists of 40 CCDs, arranged in 4 across-scan strips 
and 10 along-scan rows (the Gaia naming convention for the along- 
scan direction is rows of CCDs and columns of CCD pixels, and 
across-scan we speak of strips of CCDs and lines of pixels). In prin¬ 
ciple, the filters can be distributed among the 40 CCDs in any com¬ 
bination that is desirable from a scientific point of view. A prefer¬ 
able situation, for example for variable-star science, would be to 
have identical filters on each of the 10 CCD rows within a single 
CCD strip: this would yield quasi-simultaneous measurements in 
the different passbands on each field-of-view crossing, independent 
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Figure 2. The G and GS Gaia broad passbands corresponding to the white 
light observations in the Astro and Spectro focal planes, respectively. 

of an object’s across-scan coordinate. Such a solution would, how¬ 
ever, limit the number of BBP passbands to the number of CCD 
strips, i.e., 4. 

2.5 Spectro/MBP 

The MBP consists of 40 CCDs arranged in 20 strips and 2 rows. 
The CCDs in the first 10 strips are illuminated directly from the 
Spectro telescope. These CCDs are red-enhanced, which means 
they are thicker than ‘normal CCDs’ and have a red-optimized anti¬ 
reflection coating. The CCDs in the last 10 strips receive only the 
blue light from the Spectro telescope. These detectors are blue- 
enhanced, which means they have the same thickness as ‘normal 
CCDs’ but a blue-optimized anti-reflection coating. Four of the red- 
enhanced strips act as sky mappers and do not have filters; the as¬ 
sociated broad passband is called GS. The remaining 6 strips with 
red- and the 10 strips with blue-enhanced CCDs may be equipped 
with filters, defining up to 16 different passbands, assuming the 
two CCD rows in each strip have identical filters (for the sake of 
redundancy). One of the red passbands will cover the same spectral 
region as the RVS, around 860 nm, leaving 15 free strips. 


3 THE G AND GS PASSBANDS 

The G and GS passbands corresponding to the white light observa¬ 
tions in the Astro and Spectro instruments, respectively, are shown 
in Fig. [2] They cover the wavelength range from 400 to 1000 nm 
and 350 to 1025 nm, with the maximum energy transmission at 
~ 715 and ~ 765 nm and the full width at half maximum of 408 
and 456 nm, respectively for G and GS. The relation between the 
associated magnitudes and the Johnson V magnitude is shown in 
Fig-0 Very red objects (either intrinsically red or highly reddened) 
are much brighter in G and GS than in V, and therefore, the Gaia 
limiting magnitude of Gn m ~ 20 translates into Vji m ~ 20-25, 
depending on the colour of the observed object. 

The estimated precisions for the G-magnitudes per focal plane 
transit and at the end of the mission, computed as described in Sec- 
tion l6.2l are shown in Fig. [4] Taking into account the photon noise 
from the source, the background and the read-out noise, precisions 
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V-I c 

Figure 3. The relation between G — V and GS — V and V — Ic for 
the white light passbands in the Astro and Spectro focal planes. Every star 
in the spectral library of lPicklesj fl998 ) is represented by a filled symbol. 
Open symbols correspond to the same stars reddened by Ay= 5 mag. 
Reddening vectors run parallel to the colour-colour relationship. Pickles’ 
library was chosen because it extends to very red objects. No appreciable 
differences for stars with T e fj ^ 4000 K are obtained when using spectral 
energy distributions from the BaSeL2.2 library or other libraries. 

of ~ 10 and ~ 1 millimag are achievable at V ~ 19. This implies 
that the precision of the G measurements is ultimately limited by 
the calibration errors (see Section l 6 /il . 

Gaia will provide distances at the 10% accuracy level for some 
100-200 million stars, which combined with estimates of the G 
magnitude and the interstellar extinction, will yield unprecedented 
absolute magnitudes, in both accuracy and number. 

The G passband also yields the best signal to noise ratio 
for variability detection among all the Gaia passbands. Gaia will 
monitor millions of variable stars ( eclipsi ng binaries, Cepheids, 
RR L yrae, Mira-LPVs, etc). lEved j2005) and lEver & Mignard 
(2005) provide comparisons with other variability surveys and a 
detailed discussion of the effects of the variable time sampling and 
number of observations due to the scanning law. 


4 DESIGNING THE PHOTOMETRIC SYSTEM 

Many ground-based photometric systems (PSs) exist but none sat¬ 
isfies all the requirements of a space-based mission such as Gaia : 
portions of the spectrum limited on-ground by telluric O 3 opacity 
in the blue and O 2 and H 2 O absorption bands in the red are acces¬ 
sible to Gaia; existing PSs have usually been designed for specific 
spectral type intervals or specific objects, while Gaia photometry 
must deal with the entire Hertzsprung-Russell diagram, must al¬ 
low taxonomy classification of solar-system objects and must be 
able to identify quasars and galaxies. In addition, Gaia allows the 
extension of stellar photometry to Galactic areas where the clas¬ 
sical classification schemes may be no longer fully valid because 
of systematic variations in element abundances in stellar atmo¬ 
spheres and in interstellar matter. Finally, Gaia has to astrophys- 
ically characterize objects over a very large range of brightness, 
from G ~ 6 to the faint limit of Gn m ~ 20, and consequently the 
width of the passbands will reflect a trade-off between sensitivity 
to physical parameters and the possibility to measure faint stars. 
Therefore, designing a new PS is the best approach to ensuring that 



13 15 17 19 21 23 
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Figure 4. Estimated precision of the G magnitude per focal plane transit 
and at the end of the mission according to equation 0 assuming £T ca i = 
0. In this paper, the contribution of calibration error to the end-of-ntission 
precision is assumed to be 3 millimag (horizontal line). 


Gaia's ambitious goals will be achieved. This became clear early 
on when initial efforts to use existing PSs showed that these failed 


to cover all Gaia requirements. A fairly 

:omplete census of exist- 

ing PSs can be found in 

Straizys 

1993); Moro & Munari (2000); 

Fiorucci & Munari 1 2003 

); Bessell 

12005 

). 


Criteria for the design of a new PS have to be established a 
priori. As it is widely known, narrow passbands are very efficient 
in measuring specific spectral features, but have low performance 
for faint objects; broad passbands yield low photon-noise for faint 
objects but cannot give one-to-one determination of astrophysical 
parameters. The ultraviolet contains important information on as¬ 
trophysical parameters, but Gaia will not be able to measure faint 
red objects in this spectral range. Therefore, a compromise between 
the different options (number, location, and width of the passbands 
and their total exposure time) is needed to achieve a PS that is maxi¬ 
mally capable of separating objects with different astrophysical pa¬ 
rameters in filter flux space. Such a PS should allow both accurate 
discrete source classification (i.e. identification of stars, galaxies, 
QSOs, solar-system objects, etc.) and continuous parametrization 
(values of astrophysical parameters, e.g. for stars T e g, [M/H], etc) 
for all Gaia targets across the astrophysical parameter space. 


4.1 Principles for the design 

Over the years of mission study and development, the instrumental 
concept of Gaia has evolved. As a consequence, the requirements 
and constraints for designing the PS have evolved in terms of the 
number of passbands, the exposure time available per passband, 
the wavelength coverage in BBP and MBP, the spatial resolution of 
Astro and Spectro, the goals of BBP and MBP, the complementarity 
of BBP, MBP and RVS measurements, etc. 

With the consolidation of the study-phase of the Gaia instru¬ 
ment design, the requirements for BBP and MBP are stabilized. 
The PS accepted as the baseline for the mission and presented in 
this paper is based on the following constraints: 

• The instrument design is as described in Section[2] 
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• Photometry from Astro has to account for chromaticity in or¬ 
der to achieve microarcsec astrometry, which implies a measure¬ 
ment of the spectral energy distribution of each object with con¬ 
tiguous broad passbands covering the whole wavelength range of 
the G passband (see Section l5.1.2L 

• The photometry from Spectro has to provide the astrophysical 
characterization of the observed objects. 

• The photometry from Astro has to provide the astrophysical 
characterization of the observed objects in dense fields (at stellar 
densities larger than 200 000-400 000 stars/sq deg; mainly in the 
bulge, some Galactic disk areas and globular clusters). 

• One of the medium passbands in Spectro has to measure the 
flux in the wavelength range covered by RVS (848-874 nm). 

• The trigonometric parallax is available and has to be used to 
constrain the luminosity (or the absolute magnitude). 

• The information from RVS is not used for designing the PS 
because of the different limiting magnitude and spatial resolution 
of that instrument (in the actual Gaia data processing it is foreseen 
that astrometry, photometry and RVS measurements will be used 
together to derive astrophysical parameters, when possible). 

• The Gaia PS is optimized for single stars, where priority is 
given to those types crucial for achieving the Gaia core science 
case, namely the unravelling of the structure and the formation his¬ 
tory of the Milky Way. These stars are named ‘scientific photomet¬ 
ric targets’ or simply ‘scientific targets’ (ST; see Section l4~3l and 
they constitute our test population. 

• Every ST is characterized by its astrophysical parameters 
(APs), where we consider T e ff, log<?, chemical composition and 
Ay 1 . The chemical composition is described by the metallicity, 
[M/H], and the a-element enhancements, [a/Fe]. 

• The error goals for the AP determination are established for 
every ST (o+g oa i; see Section l4~2l and Appendix lAi. 

• The actual performance of a given PS with respect to the error 
goals is measured using an objective ‘Figure of Merit’ (see Sec- 
tion l4.2l and Appendix I aI. 

• The global degeneracies of the PS have to be evaluated (see 
Section EOl . 

• Additional merits such as, for example, the performance with 
respect to discrete object classification, the performance for non-ST 
objects, etc. have to be considered. 

The procedure to come to a baseline PS iBrown et alJ2004ih is 
based on the maximization of this ‘Figure of Merit’, a minimization 
of the global degeneracies and an evaluation of additional merits. 

4.2 Figure of Merit 

For a given photometric system the ‘Figure of Merit’ (FoM) is con¬ 
structed by calculating for each ST and each of its APs, pk. the ratio 
a fc, P ost/crfc,goai• Here dk ,post is the estimate of the error than can 
be achieved with the given PS forp^ and ak, g oa i is the above men¬ 
tioned error goal. The procedure for estim ating (jfc. DOS t and the defi¬ 
nition of the FoM were proposed bv lLindegrenl 12003b ). For a given 
AP, pk, and a PS with measured fluxes cf>j (j = 1 ,n passbands), 
a fc, P ost is estimated using the sensitivity of the PS to that parameter 
(d(pj/dpk ) and the errors of the photometric observations ( tj ). The 
latter are based on the noise model for the instrument-PS combina¬ 
tion (see Section l6~2l . The details are given in Appendix IaI 

The global FoM as given in equation <A6l in Appendix IaI is 

1 Throughout the paper Ay means A\ = 550 i.e. the monochromatic inter¬ 
stellar extinction at A = 550 nm 


a weighted sum of the individual FoMs of every ST, which in turn 
are weighted sums of the ratios o+ pos t/<7fc,g oa i for each of the APs 
Pk (equation IA4> . A higher value of this FoM indicates a better 
performance for the given PS. The global FoM thus describes the 
performance of a PS across the HR-diagram by taking into account 
the errors that can be achieved and the relative priorities of the dif¬ 
ferent STs and their APs. The local degeneracies in the AP deter¬ 
minations (i.e., correlations between the errors for different APs for 
a given ST, see Section l4~5l are also taken into account in the FoM 
formalism. Each local degeneracy between APs leads to increased 
standard errors crfc, pos t and thus a lower FoM. 

In our implementation the derivatives d<f>j/dpk and the error 
estimates <Tfc iPOSt are calculated numerically from simulated pho¬ 
tometric data (Section The calculation thereof requires (syn¬ 
thetic) spectral energy distributions (SEDs) of the STs and a noise 
model for the photometric instruments. The matrix B from equa¬ 
tion <A3> includes the a priori information for the AP vector p and 
in our case corresponds simply to the range of possible values of 
the parameters pk . The information from the parallax and its error 
is incorporated into B following fLindegreS 1200 4). Other informa¬ 
tion could be added (known reddening in a certain Galactic loca¬ 
tions, ranges of abundances according to Galactic population, etc) 
but this will introduce our preconceptions of the structure and stel¬ 
lar populations of the Galaxy into the PS design and we therefore 
did not include such constraints. 

The MBP data will have a lower spatial resolution than the 
BBP data. For non-crowded regions on the sky (with respect to 
MBP) we assume that the MBP and BBP photometry will always 
be combined for the estimation of APs. Hence the calculation of 
the achievable posterior errors is always done by combining BBP 
and MBP data. For dense stellar fields, only BBP will be available 
and in that case the achievable AP errors and the FoM have been 
computed using only broad-band photometry. 

Spectral energy distributions of the STs_in the test popula- 
tion were taken from th e BaSeL2.2 ^Leieum^yil^J9^), NextGen 
jllauschildt et alJlll999t) and AIARCS etal]f2003ib li¬ 

braries. The BaSeL2.2 librar y is a compilati o n of synthetic spectra 
from libraries p ublished by [Kutu ctI (l979t), iBessell et alJ 1 1 9891) . 
iFluksU 19941) and lHauschildt et alJll999li ~lt covers the whole HR- 
diagram and [M/H] abundances from —5 to +1 dex, but with 
solar a-element abundances. A new version of the NextGen li¬ 
brary (NextGen2) has been built taking into account Gaia mission 
needs faauschildt et all2003l ). This library includes SEDs for stars 
cooler than 10 000 K with [M/H] ranging from —2 to +0 dex and 
[a/Fe] fr om —0.2 to +0.8 dex. A more extended grid is currently 
available iBrott &Hauschildtf2005 ). Finally, a new version of the 
MARCS library taking into account non-solar a-element abun¬ 
dances has been created specifically for Gaia studies. This library 
provides coverage between 3000 K and 5000 K, [M/H] from —4 
to +0.5 dex and [a/Fe] from +0 to +0.4 dex. 

Empirical libraries, like those by IPicklesI 61998/ and 
iGunn & Strvkeil 1 1 9831) . do not provide full coverage of the HR- 
diagram and the corresponding chemical composition range, and 
therefore they are not appropriate for computing the FoM values. 

Finally we would like to make a few remarks here about the in¬ 
terpretation of the calculated values of <Tfc, pos t and the correspond¬ 
ing FoM for a set of proposed PSs. The values of dk , pos t calculated 
as described in Appendix|3should not be taken as the actual errors 
that will be achieved by the Gaia photometric system. They repre- 


2 see also http://marcs.astro.uu.se 
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Table 2. The main science goals for Ga ia concerning Galactic structure and evolution studies and the main corresponding tracers, from the 'Gaia Concept and 
Technology Study Report’ IESAl2000l) . 


Science goal Main tracers 


Chemical abundance galactocentric distribution 

Galactocentric age gradients, star formation rate 

Disentangling age-metallicity degeneracy 

Star formation history 

Detailed knowledge of luminosity function 

Halo streams, age and chemical abundance determinations 

Outer halo (R > 20 kpc), accretion and merging 

Earliest phases of evolution of the Galaxy 

Distance scale 

Thick disk formation mechanism (pre- vs. post-thin disk) 
Merging and/or diffusion 

‘In situ’ gravitational potential, K z , age-velocity relation 
Large scale structure (warp, asymmetry) 

Interstellar medium distribution 

Large scale structure (spiral arms, star-formation regions) 
Star formation rate in the bulge 
Shape of the bulge, orientation, bar 


G and K giants 

HB stars, early-AGB, main sequence turn-off stars 

main sequence turn-off, sub-giants 

main sequence stars earlier than G5 V, sub-giants 

low main sequence stars 

G and K giants 

G and K giants and HB stars 

most metal-poor stars, C- sub-giants and dwarfs near the turn-off 

RR-Lyrae, Cepheids 

G and K giants, HB stars 

high velocity A-type stars 

F-G-K dwarfs 

K-M giants 

O-F dwarfs 

OB stars, supergiants 

red giant branch, asymptotic giant branch 

red giant branch, asymptotic giant branch, Red clump stars 


sent the achievable precision if the synthetic spectra represent the 
true stars and if the noise model is correct. This issue is discussed 
more thoroughly in Sectional What makes the FoM a powerful tool 
is that it enables an objective comparison of different PS proposals, 
based on a set of agreed error goals and scientific priorities. In ad¬ 
dition for each PS a detailed study can be made of its strengths and 
weaknesses as compared to other PSs by examining the FoM for in¬ 
dividual STs and for groups of STs (such as specific types of stars, 
populations in certain Galactic directions, bright versus faint stars, 
reddened versus unreddened, etc). Once the best PS proposal has 
been chosen it can be further tuned by using the FoM procedure to 
improve its sensitivity to certain APs or to improve the performance 
for certain groups of stars. This objective figure of merit approach 
has not been used before in the design of a photometric system. 


4.3 Test population and error goals 

According to the scientific goals of Gaia (see Tabled, for every 
Galactic stellar population several kinds of stars were selected as 
STs and a priority, expressed as a numerical weight, was assigned 
to them IJordi et all2004dB) . These stars have been considered in 
different directions in the Milky Way (toward the centre, the anti¬ 
centre and perpendicular to the Galactic plane) and at different dis¬ 
tances (0.5, 1, 2, 5 , 10 and 30 kpc) in accordance with a Galaxy 
model iTorra et alJI 19991) . The STs have been reddened according 
to a 3D interstellar extinction model IPrimmel et all2003l) . assum¬ 
ing a standard extinction law (Ry = 3.1). Bulge stars in areas of 
high and low interstellar extinction were also included. This results 
in a set of about 6500 targets each with a parallax er ror estimated 
from its colour and apparent magnitude lESAll2000t) . The targets 
in the halo are considered to have [M/H] abundances from —4 to 
— 1 dex and [a/Fe] abundances from +0.2 to +0.4 dex; for the 
thick disk targets [M/H] ranges from —2 to 0 dex and [a/Fe] from 
0.0 to +0.2 dex; for the thin disk targets [M/H] ranges from —1 
to +0.5 dex and [a/Fe] from —0.2 to 0.0 dex; and finally for the 
bulge targets we assume the same range for [M/H] as for the thin 
disk and [a/Fe] ranges from 0 to 0.4 dex. The ranges of tempera¬ 
tures and absolute luminosities are assigned according to the spe¬ 


cific range of ages of each stellar population. The values of T e fj 
range from 3000 to 40 000 K and the values of log g are in the 
range 0.0-5.0 dex. 

We adopted the following precision goals (0+ goa i): 

• T e ff for A-M stars: OT eff /T e s = 1 — 2% 

• T e ff for O-B stars: crT eff /Teff = 2 — 5% 

• Ay: cta v = 0.1 mag at Ay < 3.0 mag, <ja v = 0.5 mag at 
Ay > 3.0 mag 

• My for stars with cr n /n ^ 10%: assumed known 

• loggr for stars with ov/tt > 10%: criog S = 0.2 dex 

• [M/H] (not to be determined for OB stars and supergiants): 
oJm/h] =0.1 dex 

• [a/Fe]: o-[ a / Fe ] < 0.3 dex 

• AGB: Carbon/Oxygen classification 

4.4 Photometric system proposals 

In this paper we discuss in detail only the baseline PS for Gaia that 
results from the optimization process described above. However, 
many proposals for a PS for Gaia have been studied. 

Different approac hes have been used for th e definition of 
the pa ssbands.|MunariMl999|),|Gremrn_eLal . 1199^i. [straiS v s et al] 
i200d), I.Ionhetjfi200^]2004 nVansevi£iu^^^ridiiug]200T 
l2003t> and iLindeerenl 12003 J) based their BBP and/or MBP pro¬ 
posals on their astrophysical exper tise and/or the per f orman ce of 
the existing photometric systems . ITautvaisiene et alJ 120031) and 
ITautvaisiene & Ed vardssorl l2005t) provided guidelines for design¬ 
ing a PS sensitive to C,N,0 and q-p r ocess elements. W e note 
that the sets of filters bv lMunar i 1 1999 + lGrenon et alJ 1 19991) . and 
IStraizvs et all J2000) were designed with substantially different 
constraint s imp osed in early phases of the Gaia instrument design. 

[Bailer-Jones! <12004) developed a novel method for designing 
PSs via a direct numerical optimization. By considering a filter sys¬ 
tem as a set of free parameters (central wavelengths, FWHM etc.), 
it may be designed by optimizing a figure of merit with respect to 
these parameters. This figure of merit is a measure of how well the 
filter system separates the stars in the data space and of how much 
it avoids degeneracies in the AP determinations. The resulting fil- 
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ter systems tend to have rather broad and overlapping passbands 
and large gaps in between at the same time. Some commonalities 
with conventional PSs may be recognized. Although first analyses 
showed these systems to yield FoM values only slightly inferior to 
conventional PS proposals, the method was considered too novel 
and further studies were not pursued. 

iLindegreri 1200 li) proposed several BBP systems, differing in 
the number of passbands and their profiles, and tested their perfor¬ 
mance with respect to chromatic effects estimation. He concluded 
that five passbands do not provide a clear advantage over four pass- 
bands, an d that th e overlapping of the passbands is more important. 

iHeiter et~ . ( 12005 ) used the Principal Components Analysis 
technique to design a set of BBP passbands which is optimal for 
a subset of STs with effective temperatures between 3000 and 
5000 K. Photometry was simulated for a grid of sets of four broad 
filters each. The FWHMs were varied in equidistant steps, while 
the central wavelengths were determined by requiring contiguous 
filters covering the whole wavelength range. The resulting set of 
passbands was found to perform better than all other proposals for 
bulge stars but worse than all othe rs when evaluated for the com¬ 
plete set of STs et all2004ah . 

All proposed photometric systems were evaluated using the 
FoM procedure and the results were returned to the proposers 
who were then given the opportunity to refine their respective PSs 
and submit improved versions as a new PS proposal. Updates of 


IStrains et al. 2004, 

Knude & H0g 2004, H0g & Knude 2004a b 

andjjprdi & Carrasco 

2004a hi). Finally, after several trials for the 


fine tuning of a fe w individual passbands, the C1B and C1M sets 
dJordi et all2004albl) were adopted as the baseline PS for Gaia. De¬ 
tailed results of the evaluation of all PS proposals can be found in 
fjordiet'a i] f2004Jdii . 


4.5 Local and global degeneracies 

The goal of designing a photometric system that allows both ac¬ 
curate discrete classification and continuous parametrization of ob¬ 
served sources translates to demanding that the degeneracies in the 
PS are minimized. 

The discrete classification problem concerns the way very dif¬ 
ferent parts of AP space are mapped by the PS into the filter flux 
space. For example, due to extinction a highly reddened O-star may 
to first order look like a nearby unreddened cool dwarf. A good PS 
should as much as possible be free of such ‘global’ degeneracies. 
The continuous parametrization problem concerns the way a PS 
maps a small region around a certain object in AP space onto filter 
fluxes. Here it is the ‘local’ degeneracies that are important. Exam¬ 
ples include the well know degeneracy between T e g and Av and 
the difficulty of disentangling the effects of T e g and log g for a PS 
without passbands shortward of the Balmer jump. 

Our method of evaluating the proposed PSs as outlined in Sec¬ 
tion |^| focuses on the local degeneracies. Consider the gradient 
vectors which describe how the filter fluxes respond to changes 
in particular APs. For a good PS these gradients should be large 
with respect to the noise in the data and they should ideally be 
orthogonal to each other, where the orthogonality is also defined 
with respect to the noise (i.e., for the gradients with components 
1/tj x drfij/dpk , see Appendix [X}. Large gradients mean that the 
PS is sensitive to the corresponding APs while their orthogonality 
ensures that there are no local degeneracies. The FoM takes this 
into account by calculating the posterior errors on estimated APs 
using the sensitivity matrix which contains these gradient vectors 


(see Appendix [3] for more details). Small gradient vectors are re¬ 
flected in larger errors on the estimated AP. Non-orthogonal gradi¬ 
ent vectors will also lead to larger errors and to non-zero covari¬ 
ances (i.e., correlated errors) in the posterior variance-covariance 
matrix of the estimated APs. Both effects will lead to a lower FoM 
(this is explained in more detail in Appendix^}. 

The FoM calculations do not take global degeneracies into ac¬ 
count. In fact it is assumed that one already has available a good 
classification of the object to be parametrized so that the linearized 
equations from which the FoM is derived apply. The global de¬ 
generacies reflect the highly non-linear mapping from AP space to 
filter flux space and are difficult to characterize in practice. 

We attempted to compare the different PS proposals with re¬ 
spect to global degeneracies by employing self-organizing maps 
to explore how the different STs cluster in data space and how well 
they can be separated. The results were inconclusive as the different 
PS proposals showed rather similar behaviour. This is plausible be¬ 
cause to first order all the proposed filter systems were similar. The 
sets of BBP passbands were very much alike and the sets of MBP 
passbands all had a set of blue and a set of red filters with a gap 
between these sets from ~550 to ~700 nm and they all had an Ha 
filter in this gap. With a continuous sampling of stellar parameters 
a filter system will define a complex manifold in the space of filter 
flux vectors onto which each star will be mapped. It is the over¬ 
all shape of this manifold that determines the presence or absence 
of global degeneracies in the PS. Hence it may be that the filter 
systems that had been considered all define roughly the same man¬ 
ifold, the differences between filter systems only being manifest at 
the local level (where the FoM calculations are more relevant). 

Characterizing the global degeneracies will be an important 
task in the context of the automatic classification effort for the Gaia 
data processing. A full understanding of the behaviour of the Gaia 
PS is essential for setting up appropriate discrete classification and 
continuous parametrization algorithms. 


5 THE Gaia PHOTOMETRIC SYSTEM 

In this section we describe in detail the baseline photometric system 
for Gaia, called Cl, which consists of the C1B and C1M broad and 
medium passbands. The role of each of the passbands is described 
in relation to spectral features and astrophysical diagnostics. 


5.1 The C1B broad passbands 

The C1B component of the Gaia PS has five broad passbands cov¬ 
ering the wavelength range of the unfiltered light from the blue to 
the far-red (i.e. 400-1000 nm). The basic response curve of the fil¬ 
ters versus wavelength is a symmetric quasi-trapezoidal shape. The 
filters were chosen to satisfy both the astrophysical needs and the 
specific requirements for chromaticity calibration of the astromet¬ 
ric instrument (see Section 0. T he C1B set of passbands evolved 
fro m the convergence _of the IStraizvd l2004) . iLindegrer] l2003al> 
and Jordi & Carrasccl l2004al) proposals. The specifications of the 
filters are given in Table[3] Figure[5]shows the spectral response of 
the passbands. The estimated end-of-mission precisions, computed 
as described in Section lfi~2l are shown in Fig.|b| 

The Balmer discontinuity and the H/J line limit the blue and 
red edges of the C1B431 filter, respectively. The reddest filter 
C1B916 is designed to measure the light between the Paschen 
jump and the red limit of the sensitivity of CCDs in the Astro 
focal plane. The filter C1B655 is centred on the Ha line and its 
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Figure 5. Response curves of the filters in C1B folded with the optics trans¬ 
mission and the QE of the CCDs. The spectral energy distributions of solar 
metallicity AO V, G2 V and K7 V type stars in units of W m —2 Hz -1 , taken 
from the BaSeL-2.2 library, are overplotted with (black) dotted lines (the 
vertical scale is arbitrary). The dashed (blue and violet) lines correspond to 
G2 V and K7 V type stars with [M/H] = —1 (when compared with the 
dotted lines they show the effect of a change of [M/H]). The dashed (red) 
line shows the change of the spectrum of an AO type star due to a change 
of luminosity. The solid (orange) line is the N-type carbon star AW Cyg 
extracted from the lGunn & StrvkeH I19831 ) library of stellar spectra. 


width has been optimized together with C1M656 in C1M (see 
Section l5*2l . The blue and red limits of the two remaining filters 
(C1B556 and C1B768) are consequently set in order to provide full 
coverage of the wavelength range in G (i.e. avoiding gaps between 
passbands). Four passbands are enough for chromaticity calibra¬ 
tion, but five passbands are preferred to four, when the classifica¬ 
tion and astrophysical parametrization of stars is considered. Since 
only four strips of CCDs are available in BBP, the two broadest 
passbands (C1B556 and C1B768) are implemented together in one 
strip. Flence, the number of observations with these two passbands 
will be half as much as with the other C1B passbands. 

The current design of the payload foresees no ultraviolet sen¬ 
sitivity for the Astro/BBP instrument (see Section^. The near ul¬ 
traviolet is the most important for stellar classification. The Balmer 
jump is the feature in the spectra of B-A-F type stars most sen¬ 
sitive to the temperature and gravity. It also contains information 
on the metallicity of F-G-K type stars. The absence of an ultravio¬ 
let passband in BBP is compensated with the inclusion of a broad 
UV passband (C1M326) in MBP (see Section im . The classifica¬ 
tion and parametrization of objects in Gaia is done using the BBP 
and MBP measurements together and therefore the lack of a UV 
passband in BBP should not be a drawback. However, since the 
Spectro instrument has a lower angular resolution than Astro, the 
combination of BBP and MBP data is not always possible. BBP 
will be the only tool for classification of stars in the crowded fields 
with stellar densities larger than 200 000-400 000 stars/sq deg (see 
Section l63l . Such stellar densities are found in some areas of the 
bulge and of the disk iRobin et all2005l : lDrimmel et all2005l) and 
most of these areas have low interstellar extinction (such as Baade’s 
window). In dense areas the trigonometric parallax and the Paschen 
jump will provide luminosity parametrization. 


Table 3. Specifications of the filters of C1B implemented in Astro. 


Band 

C1B431 

C1B556 

C1B655 

C1B768 

C1B916 

Abiue (nm) 

380 

492 

620 

690 

866 

A re d (nm) 

482 

620 

690 

846 

966 

A c (nm) 

431 

556 

655 

768 

916 

AA (nm) 

102 

128 

70 

156 

100 

SX (nm) 

10,40 

10,10 

10,10 

10,40 

40,10 

e (nm) 

2,2 

2,2 

2,2 

2,2 

2.2 

T max (%) 

90 

90 

90 

90 

90 

^strips 

1 

0.5 

1 

0.5 

1 


^bfue. A re d - wavelengths at half-maximum transmission 
A c : central wavelength= 0.5(Abi ue + A re d) 

AA: FWHM 

<5A: edge width (blue, red) between 10 and 90% of T m ax 
e: manufacturing tolerance intervals centred on Abi ue and A re< j 
Tjnax- maximum transmission of filter 

n strips : Number of CCD strips carrying the filter: C1B556 and C1B768 
share one CCD strip 



f 4 16 18 20 22 14 16 18 20 22 
V V 

Figure 6. Estimation of the end-of-mission precisions for the five passbands 
in C1B as a function of V, computed according to equation and taking 
(T ca l = 0 mag. In this paper, the contribution of calibration errors to the 
end-of-mission precision is estimated to be 3 millimag in every passband 
(horizontal line). The errors for C1B556 are all nearly the same for a given 
V magnitude because the mean wavelengths of these passbands are very 
similar. 

5.1.1 Astrophysical diagnostics 

The response of the C1B431 filter at the shortest wavelength is 
asymmetrical, with a red edge that is less steep. This is done to 
compensate the shift of the maximum of the response function red- 
ward due to the slope of the quantum efficiency curve and the re¬ 
flectance curve for six silver surfaces. As a result, the response 
function of the C1B431 passband becomes similar to that of the 
B passband of the UBV system. The mean wavelengths of both 
passbands are also similar: 445 nm for C1B431 and 442 nm for B. 

The mean wavelength and half width of the C1B556 response 
function is very similar to that of the V passband. As a result 
the colour index C1B431-C1B556 will be easily transformable to 
Johnson’s B-V and vice versa. Analogously, the C1B768 pass- 
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band can be easily related with the Cousins /, the SDSS i' and 
the HST 814 passbands and the colour index C1B556-C1B768 is 
transformable to V—Ic, f'-% and HST 555-814. This will facilitate 
the comparison of the numerous ground-based investigations in the 
BV system and the large number of observations being done in the 
far-red passbands with the Gaia results. As can be seen from Fig.[5] 
the differences of the ‘blue minus green’, ‘blue minus red’ or ‘blue 
minus far-red’ magnitudes may serve as a measure of metallic-line 
blanketing, although with much less sensitivity than a colour index 
containing the ultraviolet. In the C1B431-C1B556 vs. C1B556- 
C1B768 diagram the deviations of F-G metal-deficient dwarfs and 
G-K metal-deficient giants from the corresponding sequences of 
solar metallicity are up to 0.07 and 0.20 mag, respectively. 

The combination of the fluxes measured in the C1B655 pass- 
band and the narrow passband C1M656 in MBP form an Ha in¬ 
dex primarily measuring the strength of the Ha line. The Ha in¬ 
dex shares the same properties as the [3 index in the Stromgren- 
Crawford photometric system. It is an indicator of luminosity for 
stars earlier than A0 and of temperature for stars later than A3, 
almost independent of interstellar extinction and chemical compo¬ 
sition. The same reddening-free index may be used for the identifi¬ 
cation of emission-line stars. 

The two remaining red and far-red passbands, C1B768 and 
C1B916, give the height of the Paschen jump which is a function 
of temperature and gravity. Although the maximum height of the 
Paschen jump is 0.3 mag only, i.e., about four times smaller than the 
Balmer jump, it still provides the needed information if its height, 
C1B768-C1B916, is measured with high accuracy (not lower than 
1%) and this is reachable up to about V ~ 17—18 (see Fig. [6}. The 
colour indices C1B556-C1B768 and C1B556-C1B916 for unred¬ 
dened late-type stars may be used as indicators for the temperature. 
These colour indices (and C1B768-C1B916) also allow for the sep¬ 
aration of cool oxygen-rich (M) and carbon-rich (N) stars. 

The C1M326-C1B431 vs. C1B431-C1B556 diagram has the 
same properties as the U — B vs. B — V di agram of John son’s 
system or similar diagrams for other systems <Straizvslll992 l. Su¬ 
pergiants are well separated from the main-sequence stars. Metal 
deficient F-G dwarfs and G-K giants exhibit ultraviolet excesses 
up to 0.4 mag. Blue horizontal branch stars show ultraviolet defi¬ 
ciencies up to 0.3 mag, while white dwarfs are situated around the 
interstellar reddening line of O-type stars. 


5.1.2 Chromaticity evaluation 

Although no refracting optics is used for the astrometric field, the 
precise centre of a stellar image is still wavelength-dependent be¬ 
cause of diffraction and its interplay with the optical aberrations 
of the instrument. Differential shifts by up to ~10% of the width 
of the diffraction image (i.e., several milliarcsec) may be caused 
by odd aberrations such as coma, even though the resolution re¬ 
mains essentially diffraction-limited. As a result, the measured cen¬ 
tres of stellar images will depend on their spectral energy distribu¬ 
tions, and a careful calibration of the effect, known as chromatic¬ 
ity, is mandatory in order to attain the astrometric accuracy goals. 
The gross spectral energy distribution of each observed target is 
therefore needed in the wavelength range of the astrometric CCDs; 
moreover, these data are needed with the same spatial resolution 
as in the astrometric field. As stated before, the BBP set of pass- 
bands was designed with this requirement in mind, as well as on 
astr ophysical gr ounds. 

iLindegren l l2003dt showed that for the chromaticity calibra¬ 
tion, near-rectangular filters are acceptable and that the choice of 
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Figure 7. An example of the colour-dependent shifts (chromaticity) of stel¬ 
lar images that may be obtained in the astrometric field as the result of a 
coma-like optical aberration (see text for details). The shift relative to the 
geometric image centre is plotted versus a V — I c - like colour index for a 
range of synthetic stellar spectra without extinction (filled circles) and for 
Ay = 2 mag (open circles). 




Redshift z 

Figure 8. The chromatic shifts of a quasar image before (dash-dotted curve) 
and after (solid) correction for the effect calibrated by means of stellar spec¬ 
tra. A standard quasar spectrum was assumed to be observed at different 
red-shifts, and subject to the same aberrations as in Fig.0 (The dash-dotted 
curve has been displaced 4000 microarcsec downwards in the diagram to 
offset the overall shift caused by the coma). 


the separation wavelengths is more important than the edge widths. 
The author concluded that the use of four broad passbands cover¬ 
ing the wavelength range of the astrometric G passband should be 
enough to match the chromaticity constraints (r.m.s. contribution to 
the parallaxes < 1 microarcsec). 

Following the design of the passbands in ClB, a more detailed 
evaluation of residual chromaticity effects was performed. A worst- 
case scenario, representing an extreme amount of coma, was con¬ 
sidered where the wavefront aberration consisted of a third-degree 
Legendre polynomial in the normalized along-scan pupil coordi¬ 
nate, with r.m.s. wavefront error (WFE) 45 nm. Polychromatic im- 
ages in G wer e generated for a library of synthetic stellar spectra 
dMunari et alJl2005h . including some reddened by interstellar ex¬ 
tinction. Image centres were computed through a modified form of 
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Tukey’s biweight formula IPress et a ill992j) . with properties sim¬ 
ilar to the maximum-likelihood estimator to be used with the real 
Gaia data. 

FigureQshows the stellar image shifts for a range of synthetic 
stellar spectra plotted versus one of the BBP colour indices. For the 
particular WFE assumed in this example, there is a general shift of 
the image centroid by ~ 4 milliarcsec caused by the coma. FIow- 
ever, it is only the variation of this shift with the spectral composi¬ 
tion that is of concern here. The r.m.s. variation of the shift versus 
colour index is 415 microarcsec. Using linear regression against 
the synthetic stellar C1B counts <f>j (j = 1... 5, normalized to 
JW (f>j = 1), the shifts could be reproduced with an r.m.s. residual 
of 7 microarcsec. Further averaging between the ~800 astrometric 
CCD observations that are combined in a single parallax will re¬ 
duce the astrometric effect of the chromatic residuals by a factor 
0.02-0.2 depending on the degree of correlation among the indi¬ 
vidual observations, thus leading to a chromatic contribution to the 
parallax errors of 0.14-1.4 microarcsec. Since these numbers are 
based on a worst-case assumption for the WFE and a somewhat 
simplistic calibration model, it is reasonable to expect that a resid¬ 
ual contribution of 1 microarcsec can be achieved, and in particular 
that the chosen C1B passbands provide sufficient information on 
the spectral energy distribution within the G passband for this pur¬ 
pose. The residual effect is thus small compared with the statistical 
errors from photon noise and other sources even for bright stars. 

The chromaticity correction based on an empirical calibration 
against stellar BBP fluxes will be less accurate for objects with 
strongly deviating SEDs. Prime examples of this are the quasars, 
which may exhibit strong emission lines at almost any wave¬ 
length depending on red-shift. Quasars are astrometrically impor¬ 
tant for establishing a non-rotating extragalactic reference system 
for proper motions; thus chromaticity correction must work also 
for these objects. This was investigated by applying the correction 
derived as described above from stellar spectra to synthetic quasar 
images, calculated from a mean quasar spectrum observed at red- 
shifts in the range z = 0 to 6. Figure[8]shows the uncorrected image 
shift as function of z (dash-dotted curve) together with the residual 
shift after correction (solid curve). The r.m.s. image shift decreases 
from 170 microarcsec before to 29 microarcsec after correction. 
The residual curve shows artifacts that are clearly attributable to 
the limited sampling of the spectral range. For example, the nega¬ 
tive slopes for 2 = 2.3-2.8 and 2 = 3.1-3.9 correspond to the red- 
shifted Lyman-a? emission line moving through the C1B431 and 
C1B556 passbands, respectively. 

With similar assumptions as for the stars concerning the sta¬ 
tistical averaging of the effect when propagating to the astrometric 
parameters, the residual effect for the quasars will be a few mi¬ 
croarcsec. This is acceptable since these are mostly faint objects 
with much larger photon-statistical errors. 

Note that the data in Figs. □ED only represent an example 
of the image shifts that may occur. The actual behaviour depends 
strongly on the shape and size of WFE, which vary considerably 
across the field of view, and on the detailed centroiding algorithm. 


5.2 The C1M medium passbands 

The C1M component of the Gaia PS consists of 14 passbands^ and 
evolve d from the conv e rgence o f the proposals by iGrenonetalJ 
1 19991). iKnude & Hdd pOOjlkJ ^nseyicii^&BridliuTT^OO^r 
IStraizYS_etaljj200j)_and|]or^&Canasc(T l20'04ti) . The guidelines 
bv 1^utvaiIiMi^^^dvarfssoin2002l) for a-elements abundance 
determination were taken into account. The basic response curve 
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Figure 9. Same as Fig.[5]for C1M. Top and bottom figures show the 'blue' 
and ‘red’ passbands, respectively. 


of the filters versus wavelength is a symmetric quasi-trapezoidal 
shape. Their parameters are listed in Tableland the response of 
the corresponding passbands is shown in Fig.^ Six strips with red- 
enhanced CCDs are available for MBP and six red passbands have 
been designed, implemented as one filter for each strip, with the 
only constraint that MBP has to measure the flux entering the RVS 
instrument (see Section HOY For the blue passbands eight filters 
are implemented on the ten strips with blue-enhanced CCDs. Two 
strips have been allocated to each of the two UV passbands to in¬ 
crease the signal-to-noise ratio (SNR) of the measurements. 

The primary purpose of the medium passbands is the classi¬ 
fication and astrophysical parametrization of the observed objects 
(stars, QSOs, galaxies, solar-system bodies, etc.). In the case of 
the stars the goal is to determine effective temperature T e fj, grav¬ 
ity loggr or luminosity My, chemical composition [M/H], [a/Fe] 
and C/O abundances, peculiarity type, the presence of emission, 
etc., in the presence of varying and unknown interstellar extinction. 
Taxonomy classification for solar-system objects and photometric 
red-shift determination for QSOs are also aimed for. 
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Table 4. Specifications of the filters in C1M implemented in Spectro. 


Band 

C1M326 

C1M379 

C1M395 

C1M410 

C1M467 

C1M506 

C1M515 

Abiue (nm) 

285 

367 

390 

400 

458 

488 

506 

A red (nm) 

367 

391 

400 

420 

478 

524 

524 

A 0 (nm) 

326 

379 

395 

410 

468 

506 

515 

A A (nm) 

82 

24 

10 

20 

20 

36 

18 

5X (nm) 

5 

5 

5 

5 

5 

5 

5 

e (nm) 

2,2 

2,2 

2,1 

1,2 

2,2 

2,2 

2,2 

7m ax (%) 

90 

90 

90 

90 

90 

90 

90 

Type of CCD 

blue 

blue 

blue 

blue 

blue 

blue 

blue 

^-strips 

2 

2 

i 

i 

i 

i 

i 

Band 

C1M549 

C1M656 

C1M716 

C1M747 

C1M825 

C1M861 

C1M965 

Abiue (nm) 

538 

652.8 

703 

731 

808 

845 

930 

A red (nm) 

560 

659.8 

729 

763 

842 

877 

1000 

A c (nm) 

549 

656.3 

717 

747 

825 

861 

965 

A A (nm) 

22 

7 

26 

32 

34 

32 

70 

<5A (nm) 

5 

2 

5 

5 

5 

5 

5 

e (nm) 

2,2 

1,1 

2,2 

2,2 

2.2 

2.2 

2.2 

Tmax (%) 

90 

90 

90 

90 

90 

90 

90 

Type of CCD 

blue 

red 

red 

red 

red 

red 

red 

^-strips 

i 

i 

i 

i 

i 

i 

i 


-^bluei A r ed ■ wavelengths at half-maximum transmission 
A c : central wavelength= 0.5(Abi ue + A re d) 

AA: FWHM 

<5A: edge width (blue, red) between 10 and 90% of T max 

e: manufacturing tolerance intervals centred on Abi ue and A re< j 

T max : maximum transmission of filter 

^strips : Number of CCD strips carrying the filter 



10 12 14 16 18 20 10 12 14 16 18 20 10 12 14 16 18 20 22 
V V V 


Figure 10. Estimation of the end-of-mission precisions as a function of V 
for several photometric indices constructed from the passbands in C1M, and 
computed according to equation |5J with cr ca i = 0 mag. In this paper, the 
contribution of calibration error to the end-of-mission precision is estimated 
to be of 3 millimag in every passband. This is indicated by the horizontal 
line (which here includes a factor of \/2). 

5.2.1 Astrophysteal diagnostics 

The filter at the shortest wavelength is C1M326 with wavelengths at 
half-maximum transmission of 285 and 367 nm (thus it is of broad 
passband type although implemented in Spectro). Below 280 nm. 


strong absorption lines, metallicity dependent, are present already 
in A type stars. The interstellar extinction increases rapidly be¬ 
low 280 nm and reaches a maximum at 218 nm. In space the UV 
passband can be extended down to 280 nm which improves the 
determination of [M/H] for F-G-K stars because of the presence 
of many atomic lines, ionized or of high excitation. For late G 
dwarfs, the line blocking in the extended UV passband is about 2.7 
times larger than in the violet 376-430 nm domain. The red side of 
C1M326 is set at the Balmer jump. The colour indices C1M326- 
C1B431 or C1M326-C1M410 give the height of the Balmer jump, 
which is a function of T e fr and log g in B-A-F stars. In F-G-K stars 
these colour indices measure metallic line-blocking in the ultravio¬ 
let which can be calibrated in terms of [M/H]. 

The C1M379 passband is placed on the wavelength range 
where the lines corresponding to the higher energy levels of the 
Balmer series crowd together in early-type stars. The integrated 
absorption in these lines is very sensitive to logg (or My). For 
late-type stars the position of this passband coincides with the 
maximum blocking of the spectrum by metallic lines. Hence the 
colour index C1M379-C1M467 is a sensitive indicator of metallic¬ 
ity. Analogs in other photometric systems are the P magnitude in 
the Vilnius system and L in the Walraven system. 

The C1M395 passband is introduced mainly to measure the 
Call H line. The index C1M395-C1M410 shows a strong corre¬ 
lation with W(CaT*), the equivalent width of the Calcium triplet 
measured by the RVS instrument, corrected for the influence of the 
Paschen lines. The C1M395-C1M410 vs . W(CaT*) may be used 
as a log g estimator iKaltcheva et al .te'ooi Knude & Carrasco, priv. 
comm.). The C1M395-C1M410 vs. W(CaT*) plane is particularly 
useful since the effect of reddening on the colour index is only mi- 
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nor due to the small separation of the two passbands. Additional 
uses of the C1M395 passband are in assisting the [a/Fe] determi¬ 
nation and in the identification of very metal-poor stars. 

The violet C1M410 passband measures the spectrum intensity 
redward of the Balmer jump. In combination with C1M326, it gives 
the height of the jump. For K-M stars it is the shortest passband 
which when combined with longer passbands can provide temper¬ 
atures and luminosities of solar metallicity stars in the presence of 
interstellar reddening (i.e. when the stars are too faint in the ultravi¬ 
olet). Its analogs are: v in the Stromgren system, B i in the Geneva 
system and X in the Vilnius system. 

The blue C1M467 and green C1M549 passbands measure do¬ 
mains where the absorption by atomic and molecular lines is mini¬ 
mal. The flux in these domains correspond to a pseudo-continuum. 
The colour indices C1M467-C1M549, C1M467-C1M747 and 
C1M549-C1M747 may be used as indicators of the temperature 
for stars of all spectral types. The analogs of C1M467 are: b in the 
Stromgren system, B 2 in the Geneva system, and Y in the Vilnius 
system. The analogs of C1M549 are: y in the Stromgren system 
and V in the Vilnius system. 

The green C1M515 passband is placed on a broad spectral de¬ 
pression seen in the spectra of G and K type stars and formed by 
crowding of numerous metallic lines. Among them, the strongest 
features are the Mg I triplet and the MgH band. The depth of this de¬ 
pression, the intensity of which reaches a maximum around K7 V, 
is very sensitive to gravity, being deeper in dwarfs than in giants. 
The same passband is also useful for the identification of Ap stars 
of the Sr-Cr-Eu type. The same passband ( Z ) is used in the Vilnius 
photometric system. 

The C1M506 is much broader and includes the C1M515 pass- 
band region. The combination of both provides an index which 
is almost reddening free and its combination with the contiguous 
pseudo-continuum passbands (C1M467 and C1M549) provides an 
index sensitive to Mg abundances and gravity. If the luminosity is 
known from parallax. Mg abundances can be determined. The Ca II 
and Mg I spectral features _show inverse behaviour when [M/H] 
and [a/Fe] change iTautvaisiene & Edvards soi3l2002t) . and hence, 
indices using C1M395 and C1M515 allow the disentangling of Fe 
and a-process element abundances. 

The narrow passband C1M656 is placed on the Ha line. As 
mentioned in Section |5~T1 the Ha=ClB655-ClM656 index is a 
measure of the intensity of the Ha line, yielding luminosities for 
stars earlier than AO and temperatures for stars later than A3. The 
index is most useful for identification of emission-line stars (Be, 
Oe, Of, T Tau, Herbig Ae/Be, etc.). 

C1M716 coincides with one of_the de epest TiO absorption 
bands with a head at 713 nm IWahlgren et 71200.4 while C1M747 
measures a portion of the spectrum where the absorption by TiO 
bands is minimum. So, the index C1M716-C1M747 is a strong in¬ 
dicator of the presence and intensity of TiO, which depends on tem¬ 
perature and TiO abundance for late K and M type stars. For earlier 
type stars, both passbands provide measurements of the pseudo¬ 
continuum. Earlier PS proposals for MBP considered the inclusion 
of a filter centred on the TiO absorption band at 781 nm. FoM com¬ 
putations and analysis of the <7fc, poa t values showed that the AP de¬ 
termination improves by more than 10%, for [Ti/H] and T e g, if the 
passband is centred on 716 nm instead of on 781 nm. A similar but 
narrower passband has been used in the Wing eight-colour far-red 
system. 

The passband C1M825 is designed to measure either the 
continuum blueward of the Paschen jump (hence its limitation at 
842 nm for the red side mid-transmission wavelength) or the strong 


CN band for R- and N-type stars. For M stars, C1M825 measures a 
spectral domain with weak absorption by TiO. The distinction be¬ 
tween M and C stars is realized with all red passbands. At a given 
temperature, the fluxes are similar in the C1M747 and C1M861 for 
O-rich stars (the M sequence) and for C-rich stars (the C sequence), 
but very different in the C1M825 and C1M965 passbands, namely 
because of strong CN bands developing redward of 787 nm. The 
separation between M and C stars is possible even if they are heav¬ 
ily reddened. 

Similarly, C1M965 measures the continuum redward of 
Paschen jump (and in combination with C1M825 yields the height 
of the jump) or strong absorption bands for R- and N-type stars (see 
Fig.[9]bottom). Having a passband at these very red wavelengths at 
the egde of the CCD QE curve improves the interstellar extinction 
determination, which was proven through the FoM computations. 

The C1M861 passband, in between C1M825 and C1M965, is 
constrained by the wavelength range of the RVS instrument (i.e. 
848-874 nm) and hence includes the Ca IR triplet. The measure¬ 
ment of the flux of the star in this passband will help the RVS 
data reduction. The index C1M861-C1M965 measures the gravity- 
sensitive absorption of the high member lines of the Paschen series. 

Finally, the indices C1M825-C1M861 and C1M861-C1M965 
are a sensitive criterion for the separation of M-, R- and N-type stars 
(see Fig.[9]bottom). 


6 PHOTOMETRIC PRECISION 

This section describes the simulation of photometric fluxes that will 
be measured by Gaia, the associated magnitude errors and the ef¬ 
fect on these errors of crowded regions and the calibration of the 
photometric measurements. 


6.1 Simulated photometry 

For the simulation of synthetic white light ( G , GS), C1B and C1M 
fluxes and their corresponding errors, a photometry simulator Ga- 
iaPhotSim 3 was created. This tool predicts the number of photo¬ 
electrons per unit time (or per CCD crossing) for every given 
passband, observed target and Gaia instrument specifications. Ga- 
iaPhotSim also estimates the associated magnitude and the magni¬ 
tude error per transit and at the end of the mission. The inputs are: 
the instrument parameters, the SED — from an empirical or syn¬ 
thetic spectral library — of the observed object, its apparent mag¬ 
nitude and radial velocity, the extinction law and the Av value, the 
brightness of the sky background, and optionally, the coordinates 
of the object on the sky. The coordinates allow the derivation of the 
actual number of observations according to the no minal scanning 
law. The adopted interstellar extinction law is that of lCardelli et alt 
{1 9891) and the ratio of total to selective absorption can be chosen. 
The latter has a default value of 3.1. Simulations can be done for 
single and multiple stars, emission-line stars, peculiar stars, solar- 
system objects, QSOs, etc, actually for any given SED. 

To compute the object flux Sj, measured in a given photo¬ 
metric passband j and collected after a single CCD crossing, the 
following ingredients are needed: the object SED, the assumed in¬ 
terstellar extinction and the extinction law all resulting in N( A) (in 
units of photons m -2 s -1 nm -1 ), the transmission profile of the 

3 A web interface is available at http://gaia.am.ub.es/PWG/ 
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passband T)(A), the telescope transmittance T( A), the detector re¬ 
sponse Q(A) (the CCD quantum efficiency), the pupil area A and 
the single-CCD integration time r. The object flux is then given by: 

x,j 

Sj[e~] = A ■ t ■ / dA N(X) T( A) Tj( A) Q( A) (1) 

^min,j 

The flux Sj is converted into stan dard mag nitudes m,j using 
the absolute flux calibration of Vega jMesressierl 1995i) . with all its 
colours se t equal to zero. The SED of Vega has been modelled ac¬ 
cording to lRessell et all < 1994 (a Kurucz SED of T eS = 9550 K, 
log q = 3.95, [M/H] = -0.5). 

By default, the magnitude of the sky background is assumed to 
be at the level of the zodiacal light measured in Hubble Space Tele¬ 
scope (HST) observations and is set to V = 22.5 mag arcsec^ 2 . 
This assumption is very conservative for most of the sky. For exam¬ 
ple, HST measurements show that the sky background is V = 23.3 
mag arcsec -2 at high ecliptic latitudes. 

All simulations in this paper were run with the mean value for 
the total number of observations (see Table0 and with the default 
values for the extinction law and the sky background. 


6.2 Aperture photometry and associated errors 

As discussed in Section [2~3l during the observational process only 
the pixels in the area immediately surrounding the target source are 
sent to the ground in the form of a ‘window’. In most cases the pix¬ 
els in the window are binned in the across-scan direction so that the 
resulting data consists of a one dimensional set of number counts 
per sample. Following an ‘aperture photometry’ approach, it is as¬ 
sumed that the object flux Sj within a given passband j is measured 
in a rectangular ‘aperture’ of n s samples within the window. Some 
light loss can be produced due to vignetting and/or the finite extent 
of the ‘aperture’. Hence the actual measured flux will be / aper • Sj, 
where / aper ^ 1. 

The end-of-mission magnitude error (<r mj ) is computed tak¬ 
ing into account: (1) the total detection noise per sample r, which 
includes the detector read-out noise, (2) the sky background con¬ 
tribution bj assumed to be derived from rib background samples, 
(3) the contribution of the calibration error per elementary obser¬ 
vation <j ca i, and (4) the mean total number of observations n e g = 
n 0 bs x Rstri P s ■ The resulting magnitude errors are increased by a 
20% (m = 1.2) ‘safety margin’ to account for sources of error not 
considered here (such as a dependence of the calibration error on 
sky density and source brightness, see below): 

a m,j [mag] = m • }_ ■ \ a 2 a t + (2.5 • log 10 e- 

V ^eff \ ' 

[/aper ' Sj (pj T ) (l Tl 

f&per ' Sj 

The attainable precisions are shown in Figs. 0 and |6| where 
the estimated a m ,j are plotted as a function of V and spectral type 
for the G and C1B passbands, respectively. The associated end- 
of-mission errors for several colour indices formed with the C1M 
passbands are shown in Fig. 1101 A precision of 0.01 mag is obtained 
at V ~ 18 and 16 for the BBP and MBP passbands, respectively. 

For high SNR, the end-of-mission magnitude error is limited 
by a ca i. Following the approach in the Gaia Study Report IesaI 
ti 2000l) . p. 264) and evaluating the number of involved calibration 
parameters and the number of available stars for calibration pur¬ 
poses, (T ca i ranges from 0.3 to 2 millimag depending on the instru¬ 
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Figure 11. Various single transit error estimates: the dashed (blue) line is 
for the case of aperture photometry (equation the solid (green) line for 
the case of PSF fitting photometry using accurate positional information 
from astrometric processing in a non-crowded field, and the shorter (yel¬ 
low) line on the right above the green line shows the same for a crowded- 
field (400 000 stars/sq deg). The black dots are the formal errors from the 
PSF fitting solution in the crowded case, (msgo is the magnitude for a rep¬ 
resentative MBP passband at 590 nm). 



ment and the passband. Thus potentially sub-millimagnitude pho¬ 
tometric precisions can be achieved at the end of the mission. How¬ 
ever, since there is no detailed calibration model for the photomet¬ 
ric data processing at the moment, we prefer to be very conservative 
and we assume an ad hoc calibration error floor of 30 millimag per 
single observation, yielding a minimum calibration error of about 
3 millimag at the end of the mission (a m j -S 3 millimag). 


6.3 Crowded areas, image restoration 

An important limiting factor for the photometric accuracy is the 
distribution of the background flux around each source. Part of the 
background consists of diffuse emission from the general sky back¬ 
ground and zodiacal light. In addition there will be discrete back¬ 
ground sources which may be fainter than the Gaia survey limit 
(Gii m ) and PSF features of neighbouring bright stars may also con¬ 
tribute to the background. In Astro there is the added complication 
that the fields of view from two telescopes overlap on the focal 
plane and that the two overlapping fields will be different for dif¬ 
ferent scan directions. 

This background problem becomes particularly severe in the 
more crowded areas on the sky (above ~ 10 s stars/sq deg) where 
the limited resolution of the MBP instrument will lead to blend¬ 
ing of sources in addition to the effects mentioned above. In this 
case the data reduction will require a careful deconvolution and it 
was shown byi van 3 <2004 that this is possible for densities up to 
about 2-4 x 10 s stars/sq deg to magnitude 20, by using the accu¬ 
rate positional information from the astrometric processing and the 
knowledge of the PSF in Spectro. Moreover, the study also showed 
that the precisions estimated from equation Q are pessimistic (see 
Fig. ll II . In non-crowded regions, PSF photometry yields lower er¬ 
rors than those estimated with aperture photometry. 

The data processing for crowded regions can be improved 
further if we also have knowledge of disturbing sources around 
each object to a magnitude limit that goes fainter than Gn m . Im¬ 
ages providing this knowledge can in principle be reconstructed 
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from windows that were obtained for different scan directions. 
This reconst r uction process has been studied bv lNurm il EpoH and 
iDollet et all 120051) who present two different image restoration 
methods. Both studies show that by combining the Astro sky map¬ 
per windows or a longer AF window one can obtain images at a 
resolution of ~ 0.1-0.2 arcsec that go 2-3 magnitudes deeper than 
Glim- i urm I B in particular shows that using a longer AF win¬ 
dow one can map the disturbing sources in the immediate surround¬ 
ings of objects (within a ~ 2.5 arcsec diameter) to V ~ 24 for 
brightness differences A V < 8. The image reconstruction will ob¬ 
viously also be useful in identifying components of multiple stars. 

6.4 Calibration 

The photometric precision and accuracy will ultimately be limited 
by how accurately the data are calibrated. At an elementary level 
the goal of the photometric data reduction is to accurately extract 
from each CCD image the object number counts, and to transform 
these into calibrated fluxes on a standard scale. 

The effects that influence the measured flux can be broadly di¬ 
vided into four categories: (1) the mirror and filter reflectivity and 
transmission profiles; (2) the details of the CCD response, includ¬ 
ing non-linearities and charge transfer inefficiency effects; (3) the 
point spread function and uneven image motion during the transit; 
(4) the sky background and disturbance from neighbouring sources. 
The calibration process concerns the careful control and monitoring 
of these effects, all of which are time dependent. 

The Gaia observations are fully self-calibrated and need not 
rely on any extant photometric system. The variations in response 
across the focal plane, and the variations with time, will be mon¬ 
itored and corrected taking advantage of the scanning mode in 
which Gaia will be operated. For average stellar densities on the 
sky about 60 and 300 stars per CCD per second will cross the focal 
planes of Astro and Spectro, respectively. This translates to about 
660 and 3400 stars per pixel column per 6 hours (which is the spin 
period of Gaia). The same stars are observed repeatedly in different 
parts of the focal plane on time scales varying from hours to weeks 
to months, up to the mission life time of 5 years. There will thus be 
plenty of measurements to perform detailed CCD calibrations on 
scales from pixel columns to CCDs. 

The photometric system will be defined by the average re¬ 
sponse for all the CCDs over the duration of the mission. It is also 
important to establish a standard scale, where all non-linearities in 
the response are corrected. These effects will be seen as magnitude 
dependent PSFs, especially for the brighter sources. A more exten¬ 
sive discussion of the calibration issues for the photometric data 
processing can be found in tBrow ^ <2005h . 


7 PHOTOMETRY PERFORMANCES 

The role of the photometric system is to astrophysically charac¬ 
terize the observed objects, and mainly to determine APs for sin¬ 
gle stars. Thus, prior to parametrization, a classification is needed. 
Classification allows the identification of several kinds of objects 
(stars, QSOs, etc.). This is done mostly using the photometry, al¬ 
though parallaxes and proper motions will help the identification 
of extragalactic objects. Stellar parametrization with Gaia is very 
challenging due to the large range of stellar types (across the whole 
HR-diagram) encountered, the complete absence of prior informa¬ 
tion and, for the vast majority of stars, the availability of only the 19 
BBP and MBP fluxes. For some 100-200 million stars, parallaxes 


will be measured at the 10% accuracy level. For the brightest stars, 
the RVS spectra will also provide information on APs. Parallaxes 
and RVS will be used when available, although optimal exploita¬ 
tion of heterogeneous data requires sophisticated approaches. Other 
challenges are presented by: the presence of degeneracies (where, 
e.g., a set of fluxes may correspond to more than one type of star, 
see also Section EOl : the very large cosmic variance compared to 
limited data; the inevitable presence of strange objects; the prob¬ 
lem of accurately determining ‘weak’ APs ([M/H], logg) in the 
presence of dominant variance from ‘strong’ APs (T e fr, Ay). 

This section presents two different evaluations of the perfor¬ 
mance of the Gaia PS with respect to the determination of APs for 
single stars. One approach is through the analysis of the posterior 
errors introduced in Section l4~2l taking into account photometry 
and parallax information. The other approach is through the analy¬ 
sis of results obtained with parametrization algorithms specifically 
designed for the AP determination. In both cases we assume that 
the differences among the measured fluxes of the stars, and hence 
the changes in their spectra, are only due to differences of the APs 
Pk and the errors. We also assume that the SEDs, and their change 
as a function of the parameters pk, match the true stars. In other 
words, cosmic dispersion (the effects of rotational velocity, mag¬ 
netic fields, element-to-element chemical composition differences, 
differences in the interstellar extinction law, unresolved binarity, 
and so on) and inaccuracies in the SEDs libraries (absence of chro¬ 
mospheres and some opacity sources, NLTE effects, simplified ge¬ 
ometry, etc) are not taken into account. In addition, for both ap¬ 
proaches we assume that the noise model described in Section Iff2l 
is correct. Therefore, both the posterior errors as well as the errors 
estimated from the parametrization algorithms are, in principle, op¬ 
timistic. The posterior errors, which evaluate locally the sensitivity 
of the PS with respect to APs, are more optimistic than the results 
from the parametrization algorithms because global degeneracies 
are assumed not to exist (see Section 1431 . On the other hand, the 
errors on the APs estimated from the parametrization algorithms 
reflect a combination of the capabilities of the PS and of the algo¬ 
rithms themselves. Thus a large estimated error may mean a poor 
PS performance or a poor efficiency of the algorithm, or both. The 
parametrization algorithm discussed in Section l7*2l onlv considers 
the C1M passbands; information from BBP and parallax is not in¬ 
cluded. Hence, the estimated errors tend to be pessimistic. 

By the end of the Gaia mission improved stellar atmosphere 
models, improved opacity sources, improved treatment of convec¬ 
tive atmospheres, dusty atmospheres, geometry, NLTE effects, etc. 
will be available. This, supplemented with real data where appro¬ 
priate, will allow us to minimize one source of inaccuracies. 

7.1 Posterior error, ak, P ost, estimations 

The posterior errors on the APs, crk, P o st, for the baseline C1B and 
C1M PS have been computed as explained in Appendix[3]and are 
presented here. The targets from Section 1351 have been grouped 
by spectral type, luminosity class and Galactic direction and their 
mean (7T e ff.post/T^ff, (JAy, post, &iog g ,post and /h] ,post, are 

shown in Fig. 1121 as a function of distance from the Sun (i.e. as 
a function of apparent magnitude). The figure shows the error pre¬ 
dictions based on the BaSeL2.2 SED library and the combination 
of the white light, BBP and MBP fluxes with parallax informa¬ 
tion. The fluxes measured in the different passbands are assumed 
to vary only due to changes in T e g, Ay, logg and [M/H], Us¬ 
ing other SED libraries leads to slightly, but not significantly, dif¬ 
ferent o>, P o S t values. The largest differences occur for cool stars 
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(<4000 K) for which the differences among synthetic spectral li¬ 
braries are large. For these stars, discre pancies am ong theory and 
observations also exist. See lKncinskas et ail l! 2005h for a compari¬ 
son of observed colour indices and predictions by the PHOENIX, 
MARCS and ATLAS model atmosphere codes for late-type giants. 

For a given group of spectral types and luminosity classes, 
the posterior errors increase with the distance to the Sun, because 
the stars become fainter and the measurement errors of their fluxes 
increase. In the Galactic plane the stars are also increasingly red¬ 
dened (and fainter) with distance which also leads to larger errors. 
Given a distance and a group, the posterior errors increase from the 
Galactic pole direction to the Galactic centre direction due to the 
increase of interstellar extinction. For a given distance and direc¬ 
tion, the errors are larger for dwarfs than for giants and supergiants 
due to the different absolute magnitudes. 

The lower limit on the crf^post values reflects the minimum 
photometric error due to the contribution of calibration errors. As 
discussed in Section l6/il we have conservatively assumed an ad hoc 
error floor of ~3 millimag at the end of the mission. 

Bulge stars are generally faint because of their distance and 
the high extinction toward them. However, there are areas of low 
extinction, such as Baade’s window, but there the stellar density is 
very high and MBP measurements will not be available. In those 
cases the derived stellar parameters, based only on G , BBP pho¬ 
tometry and parallaxes, are less precise. Nevertheless, the estimated 
<?A V , post ranges from 0.2 to 0.6 mag, (JT eCf /T e ff from 1 to 
10% and cri oggiPOS t from 0.1 to 0.5 dex, depending on the spectral 
type. In the high extinction areas, the lower stellar densities allow 
for MBP measurements and even though the stars are very faint, the 
estimated precision of the APs is still acceptable: a a v , P ost ranges 
from 0.4 to 1.0 mag, CTT eff , P ost/7eff from 2 to 15% and eriogg.post 
from 0.3 to 0.7 dex. 

Figure fOl shows the crfc iPOS t values as a function of apparent G 
magnitude and [M/H] for F-type stars representative of isochrone 
turn-offs for the old-disk, thick disk and halo stars. This example 
has been chosen because F-type stars at the turn-off are key targets 
for many Galactic topics (see Tabled- The performance of the PS 
is excellent up to about G = 18, even for the most metal-poor stars, 
and decreases for faint magnitudes. The determination of [a/Fe] is 
only possible up to G ~ 16, while at G ss 18, the measurement 
errors are large enough to make the changes in SEDs indistinguish¬ 
able which means that <T[ a / Fe i poat ~ ofct/Fe] .prior = 0.3 dex. 

As discussed in the introduction of this Section, the 
c/c, post values are optimistic by construction and from the figures 
above they may seem too optimistic compared to the achievable 
precisions with the currently available PSs. However, they are not 
necessarily overly optimistic. Gaia photometry will provide mea¬ 
surements in 5 broad and 14 medium passbands and the G white 
light passband, to which the information contained in the parallaxes 
can be added. There is no precedent for deriving stellar astrophysi- 
cal parameters in this way. 

The at, post here represent the maximum precision that can be 
obtained. We estimate that the main uncertainty on the actual values 
comes from the imperfect knowledge of the SEDs of the true stars 
and not from the design of the photometric system itself. 

7.2 Estimating the astrophysical parameters 

The predicted errors discussed above are evaluated locally and 
hence they are based on the assumption that global degeneracies 
(e.g., the confusion between a reddened hot star and an unreddened 
cool dwarf) are already accounted for before the APs are estimated. 


(I,b) = (0 deg, 90 deg) 



Figure 12. Predicted (posterior) errors for the main astrophysical parame¬ 
ters in three Galactic directions based on G, C1B. C1M and parallax data. 
The following assumption holds for Ay: in the Galactic pole direction it 
reaches 0.3 mag at 1 kpc and is constant from this distance on; toward the 
Galactic centre it ranges from 0.3 mag at 500 pc to 10 mag at 10 kpc; toward 
the anti-centre it varies from 0.3 mag at 500 pc to 3.5 mag at 5 kpc and is 
constant for large distances. The figure legends indicate the meaning of the 
lines, where 'Class I’ means luminosity class I and ‘HB' means horizontal 
branch stars. 
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Figure 13. Predicted (posterior) errors for the astrophysical parameters as 
a function of apparent G magnitude and [M/H] for F-type subgiants. The 
top four panels show the results based on the BaSeL2.2 SED library. In the 
top right panel the metallicity increases from [M/H] = —4 at the top to 
[M/H] = +0.5 at the bottom for each vertical set of dots. The bottom left 
panel shows the predicted errors for [a/Fe] based on the NextGen2 SED 
library. The bottom right panel shows the error on [M/H] as a function of 
[M/H] for stars brighter than G = 15 (based on BaSeL2.2 SEDs). 


We now discuss the performance of the Gaia PS in a more realistic 
setting where the APs of each source have to be derived from the 
measured photometry in the presence of global degeneracies. 


7.2.1 Methods for astrophysical parametrization 

Numerous statistical methods are available for estimating parame¬ 
ters from multidimensional data. The most appropriate are super¬ 
vised pattern recognition methods, which involve a comparison of 
the observed ‘spectrum’ (the set of observed filter fluxes) with a set 
of template spectra with known APs (the ‘training grid’). 

Minimum distance methods (MDMs) make an explicit com¬ 
parison with every spectrum in the training grid to find the ‘closest’ 
match according to some distance metric (e.g. iKatz et a 0(1998). 
Because only discrete APs of templates in the training grid can be 
assigned to the observed spectrum, the precision is limited by the 
‘AP resolution’ in the grid. This may be partially alleviated by av¬ 
eraging over the APs of the nearest ne_ighbours or by interpolating 
between them (e.g. lMalvutcl2005l : lBridzius & Vanseviciusl20oH) . 

Instead of making an explicit comparison with every tem¬ 
plate spectrum, we can model the mapping between the observed 
data space, V , and the AP space, S, with a regression function 
p = /(</>; in), i.e. the APs p are a function of the fluxes <fi and 
a set of free parameters w. We then solve for the free parameters 



4.6 4.4 4.2 4 3.8 3.6 3.4 4.6 4.4 4.2 4 3.8 3.6 3.4 

log(T eff ) log(T ell ) 

Figure 14. The distribution of the astrophysical parameters in the training 
(right) and test set (left). In total there are 61 941 stars (AP combinations) 
in the training set and 9229 different stars (each represented by ten different 
noisy filter fluxes for a given magnitude) in the test set. 


w via a numerical minimization of the errors in the predicted APs 
of the training data. As both V and S are multidimensional, and 
because V —> 5 is an inverse mapping (and hence may show 

degeneracies), the problem is complex iB ailer-.Tone sl2003ll2005l) . 
Nonetheless, it offers advantages over the direct template matching, 
such as less dependence on a dense template grid, provides contin¬ 
uous AP estimates and much faster application times to new data 
(because the mapping function is learned in advance). There are 
many ways to implement a multidimensional regression, including 
neural networks (NNs) and adaptive splines (see e.g. lHastie et alJ 
120011 ) . 

Various methods have been tested by the Gaia classification 
and photometry working groups and each have their advantages and 
disadvantages. For more details see Brown 1 20031 ). In the rest of this 
Section we focus on one particular method. It must be stressed that 
this is just a preliminary investigation into AP estimation with Gaia 
and is restricted to just the C1M passbands. 


7.2.2 A regression model 

We applied a feed-forward NN to a set of simulated photometry 
in the medium-band C1M component of the Gaia PS. The train¬ 
ing grid consists of 61941 simulated ‘spectra’ showing variance in 
the four APs T eS (2000 K to 50 000 K), [M/H] (-5.0 to +1.0), 
log g (—1.0 to +5.5), Ay, the line-of-sight interstellar extinction 
(0.0 to +5.0). The source spectra were taken from the BaSeL2.2 li¬ 
brary. The distribution over the APs is discretized and non-uniform 
(Fig. 1 1 41 right panel). We used the STATNET NN code iBailer-Jones! 
1998, 2000) with two hidden layers each comprising 25 hidden 
nodes. STATNET is trained with a conjugate gradient optimizer with 
weight decay regularization to avoid over-fitting. We assess the per¬ 
formance by applying the network to a separate set of data; the 
distribution of the APs in this ‘test set’ is shown in Fig. 1141 (left 
panel). In total there are 9229 different stars (AP combinations) in 
the test set, each represented by ten different noisy passband fluxes 
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Figure 15. Parametrization errors (as defined in equation 0 for the four 
astrophysical parameters (APs) as a function of T e ff for different subsets 
of Ay and [M/H]. For a given T e ff a data point lies at the average of 
a representative temperature interval. Other than for the log g error plot, 
results are shown averaged over the range log g E [3.5, 5.5] dex, although 
the results are not significantly different for giants. For Ay and T e ff the 
results are shown for [M/H]E [—1,0] dex, for metallicity and gravity, Ay 
is limited to [0, 0.5] mag. The errors for T e ff are the fractional errors. While 
we show results for these limited AP ranges for clarity, it is important to 
realize that the model was trained and tested on the full range of APs shown 
in Fig. 114-1 there was no prior restriction on the APs of a star presented to 
the model, even though several combinations of G and APs may not exist in 
the Galaxy. All results for the C1M passbands for stars at G=15 mag using 
simulated photon fluxes are representative of the end-of-mission. 


for a given magnitude. Separate models were trained and tested at 
magnitudes G=15, 18 and 20, i.e. with the appropriate amount of 
simulated noise, from which we can assess performance as a func¬ 
tion of magnitude. All combinations of individual values of the four 
APs and the apparent G magnitude have been used without regard 
to the reality of such combinations. Moreover, we do not build in 
any prior information concerning the probability of occurrence of 
the various AP combinations. 


7.2.3 Results 

The results are summarized in terms of the average absolute errors 
1 N 

E = — ■ ^ | computed(p) - true(p)| , (3) 

p=i 


although the results are not significantly different for giants. A 
complete o vervi ew of the results for all stellar types is given in 
Will emsen et alJ 12005 all . 

For hotter stars (A, B type stars) the errors in [M/H] and T e s 
are larger than those for cooler stars (G and K type) while the Ay 
and log g errors are smaller. This behaviour is expected and reflects 
what we know of how the APs are expressed in the SED as a func¬ 
tion of the APs (in particular T e g). For stars of G type and earlier, 
the extinction is determined quite well. For example for A stars we 
find errors of 0.03 to 0.06 mag for G=15 and 0.07 to 0.15 mag for 
G=18 mag. For cooler stars (K and M), the performance degrades 
slightly, especially for highly extinct cool stars. 

Metallicity determinations are possible for cool and intermedi¬ 
ate temperature objects (M to F type) with precisions ranging from 
0.1 to 0.4 dex, even at low metallicities. For the same stars, [M/H] 
can be determined to 0.1-0.2 dex down to —2 dex metallicity. Here, 
only stars with T e ff^3500 K have been included. For stars below 
~4000 K, the estimated uncertainties may be unrealistic due to the 
inaccuracies of the current atmosphere models, which in this work 
are assumed to match the reality, as discussed previously. 

The precision of the gravity determination depends on the 
temperature, but is basically independent of metallicity. We find 
errors of E(log g ) = 0.08 to 0.4 dex for these stars at G = 15 mag 
and 0.2 to ~ 1.0 dex at G = 18 mag. For hot stars (A, B type), 
we find errors of ~0.1 dex or better, even for G = 18 mag. This 
will provide good discrimination between distant halo HB stars and 
nearer A dwarfs in those cases where the parallaxes are poor. 

Temperatures can be determined to between 1% and 5% for 
M to A type stars (G = 15 mag) and to precisions of 2% to 11% 
at G = 18 mag. There is some dependence on the extinction. For 
lower values of the extinction (AyG [0, 0.5] mag) we find temper¬ 
ature errors of E(T e ff) ^ 4% for (G = 15) and ^ 5% (G = 18) 
for all types of stars. 

We performed additional tests to determine the [a/Fe] abun¬ 
dance using the NextGen2 SED library. This increases the dimen¬ 
sionality of the AP space (from four to five) and thus its complex¬ 
ity. We used the same NN approach and the C1M passbands. We 
found errors of E([a/Fe]) ~ 0.1 dex for low to intermediate tem¬ 
perature stars at G = 15 mag. Our results also suggest that pre- 
cisions_of 0.2 dex sho uld be po ssible for G < 16.5 mag. See 
IWillemsen & Bailer-Jonesl l2005f ) for more details. 

Comparative tests were also made with MDM and Radial Ba¬ 
sis Function Neural Networks (RBFNNs) using the same data. In 
general, feed-forward NNs yielded the best results, followed by 
MDM and then RBFNN. 

7.3 Discussion 


where p denotes the p th spectrum (star) and computed(p) is the 
parametrization output provided by the network. The errors for T e s 
are given as the fractional errors. It is well known that the error in 
any AP for a given star depends significantly not only on the value 
of that AP, but also on the value of the other APs. Therefore, re¬ 
porting a single error for each AP by averaging over the full range 
of the other APs has little meaning and also depends strongly on 
the AP distribution in the test set. We therefore present errors aver¬ 
aged over selected narrow ranges of APs. This is only for reporting 
purposes. When the network is presented with a star to classify, it 
has no prior knowledge of its APs (not even that it is restricted to 
the range in the training set). 

Representative plots of the AP errors are shown in Fig. l 1 51 
for G=15. This plot is for dwarf stars, i.e. logp G [3.5, 5.5] dex, 


These preliminary results, based on the C1M passbands, demon¬ 
strate that an automated ‘bulk’ determination of astrophysical pa¬ 
rameters is possible. 

The posterior errors predicted from the sensitivity of the SEDs 
to changes in APs are lower than the errors estimated from the NN 
results because the former are immune to global degeneracies (by 
design), whereas the NN performance is degraded by its failure to 
explicitly deal with degeneracies (which produce a non-uniqueness 
in the mapping it is trying to solve). The NN performance may be 
degraded by inefficiencies in the formulation of the problem via 
multidimensional regression. Yet, the posterior errors are based on 
C1B, C1M and parallax data, while the errors with the NN are esti¬ 
mated using only C1M passbands. Taking all into account, Figs. 1121 
H3l and fT5l show an overall agreement. 
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The actual set of training data we will use for the Gaia data 
processing will be based on improved stellar atmosphere models 
supplemented with real data. However, even then the training data 
are unlikely to represent the full range of cosmic variance in the 
APs which Gaia will encounter. 

The performance predictions from the present implementation 
of the NN discussed above are pessimistic in several ways. First, 
only the C1M passbands were used. Including data from the five 
C1B passbands will improve the results (for instance, the index 
C1B655-C1M656 helps to break degeneracies among T e g, loggr 
and Ay ). Second, the use of the parallax, when available, will 
also help as it provides information on the intrinsic luminosity (for 
known apparent magnitude and estimated extinction) which will 
help solving degeneracies between log g and chemical composi¬ 
tion. Initial tests of including the parallax as an additional network 
input confirm that this improves the estimation of the gravity, es¬ 
pecially for those stellar types for which log g is otherwise poorly 
constrained (i.e. for intermediate and lo w-tempe rature stars). It also 
improves the metallicity estimates ( Iwillemsen et al]F2005bl) . This 
decrease of global degeneracies makes the estimated errors from 
NNs more similar to the predicted posterior errors. Third, for the 
brightest stars, the RVS provides high SNR spectra for all types 
of stars. This should improve the AP estimates. Finally, and most 
importantly, we know that the method we have used is limited in 
a number of ways. For example, the data to AP mapping must be 
inferred based only on the training data. Explicitly providing infor¬ 
mation on the sensitivity of passbands or passband combinations to 
APs will help. Likewise, the NNs do not yet deal with AP degener¬ 
acy. This is a particular problem at low SNR. Overall, the predicted 
PS performances from Section l7*2l are rather conservative estimates 
of what will ultimately be possible with Gaia. 


8 SCIENCE IMPLICATIONS 

As emphasized a number of times, the main scientific goal of Gaia 
is the quantitative description of the chemical and dynamical evo¬ 
lution of the Galaxy over its entire volume. This is only possible if 
physical properties of stars, through chemical abundances and ages, 
are analyzed together with kinematics and distances. To determine 
ages and abundances with an accuracy sufficient for Galactic stud¬ 
ies, temperatures (and hence extinction) and luminosities have to 
be accurately determined as well. 

In the following subsections we outline the implications of 
the Gaia photometric system performance for addressing the sci¬ 
ence case. However, an in-depth discussion is not intended. In ad¬ 
dition, we comment on the photometric performance for QSO clas¬ 
sification. As mentioned in Sections fTI and l5~l.21 a proper iden¬ 
tification of QSOs will allow the definition of a non-rotating ex- 
tragalactic reference frame which is in dispens able for astro metric 
purposes. For other objects, we refer to lKaempf et al.1 f2005 ) who 
deal with the automatic parametrization of unresolved binary stars, 
to lKolka et alJ |200^)_ for the study of performances for emission¬ 
line stars and to lCellino et alJ 120051) for a discussion about Gaia 
photometry and the parametrization of asteroids. 

8.1 Interstellar extinction and effective temperature 

The determination of the effective temperatures, absolute magni¬ 
tudes and chemical abundances requires accurate estimates of in¬ 
terstellar extinction including the interstellar extinction law. Ide¬ 
ally the extinction estimates for individual stars are based solely on 


each star’s observables, as the extinction varies on small scales in 
any field except for the closest stars. Assuming a standard extinc¬ 
tion law, Gaia's photometry allows the determination of such indi¬ 
vidual extinction measures (see Figs. ll2M13l and n*5l . though with 
significantly larger error for the late-type stars due to degeneracy 
in the Ay and T e g determinations. Meanwhile, the determination 
of T e fi for the late-type main-sequence stars remains reliable, even 
in the presence of significant extinction, to at least G = 18. Effec¬ 
tive temperatures with similar uncertainties are achievable for red 
giants in the halo where the extinction is low. 

We can identify the main contributors to the determination of 
reddening. For example, the C1M825-C1M965 vs. the C1B655- 
C1M656 (Ha index) diagram is useful for extinction measure¬ 
ments. The C1M825-C1M965 colour index is primarily meant as 
a measure of the strength of the Paschen jump, sensitive to stellar 
type, while the Ha index is nearly reddening free. For stars where 
the Ha index is a measure of the strength of the Balmer line, i.e. 
stars earlier than about G3-G5, the unreddened main sequence dis¬ 
plays a sharp locus in this colour-colour diagram. The same Ha 
index value may, however, be measured for stars on either side of 
the Balmer maximum (A1-A3 stars), but this ambiguity is lifted by 
comparison to the C1M656-C1M965 vs. Ha index diagram. Thus 
the inclusion of the Ha index allows a classical approach to deter¬ 
mining intrinsic colours from standard curves and simultaneously 
correcting the colours for the influence of metallicity and evolution 
across the main sequence. But again, intrinsic colours only result 
when the Ha index is measured precisely, while a potential draw¬ 
back of the C1M825-C1M965 index is that these two C1M pass- 
bands are among those least affected by extinction. 

If the reddening in only one colour index is determined, then 
the extinction in other passbands would have to rely on an assumed 
extinction law. Indeed, in our performance analysis (Section^ we 
have assumed _a_standard extinction law. Deviations from this stan¬ 
dard law jFitznatrickl 1999i) . due to spatial variations in the chemi¬ 
cal composition and size distributions of interstellar dust, will con¬ 
tribute to systematic errors in the individual extinction determina¬ 
tions. The breadth of Gaia’s photometric system justifies address¬ 
ing the question of whether the extinction law itself (averaged over 
the line-of-sight to a star) can be determined from the photometry. 
Since Gaia is a comprehensive survey of all point sources, the num¬ 
ber of sources is large enough that most stellar types are sufficiently 
represented even in small fields, close to the Galactic plane, where 
extinction correction is most essential. In this case, local extinction 
laws may possibly be deduced from colour difference vs. colour 
diagrams. For example, deviant extinction laws may be identified 
from Q-like indices vs. colour diagrams. 

Gaia RVS observations will provide other measures of extinc¬ 
tion. First, there is the diffuse interstellar band at 862 nm focated 
within the wavelength range of the RVS instrument ( Katz et al i 
12004) . In addition, the equivalent width of the IR calcium triplet 
gives a rough indication of the T e g ; in combination with the vir¬ 
tually reddening free C1M395-C1M410 index the calcium triplet 
can be used as a diagnostic of intrinsic colours and thus to estimate 
the extinction. Together with Gaia parallaxes, all of these measures 
of extinction (spectroscopic and photometric) will allow the con¬ 
struction of a 3D Galactic extinction map. 

Gaia measurements will also allow a statistical approach to 
the construction of a 3D extinction map. Given the parallax, main 
sequence loci may be empirically constructed from nearby, unred¬ 
dened samples provided by Gaia. Shifting these loci to a given dis¬ 
tance with a given mean extinction identifies the main sequence in 
the observed colour-magnitude diagram (including the photometric 













20 Jordi et al 


Table 5. Distances for which the relative parallax error is ~ 10%: d Q is the value of this distance for zero interstellar extinction and d a b s is the value for an 
average Galactic plane interst ellar extinction of 0.7 mag kpc -1 . V(d Q ) and V (d a b s ) are the coiTesponding apparent V magnitudes. Parallax accuracies are 
from Table 8.4 in lFSAt <2000 1 


SP 

M v 

do (pc) 

V(d 0 ) 

^abs (P c ) 

F(d ab s) 

SP 

My 

do (pc) 

V(do) 

dabs (pc) 

V (dabs) 

B1 V 

-3.2 

20 000 

13.2 

7000 

15.7 

G8III 

0.8 

9000 

15.6 

4400 

17.1 

A0 V 

0.65 

8500 

15.2 

4500 

16.8 

K3III 

0.3 

10 000 

15.3 

4800 

17.1 

A3 V 

1.5 

7000 

15.7 

3800 

17.1 

M0 III 

-0.4 

13 000 

15.2 

5500 

17.2 

A5 V 

1.95 

6500 

16.0 

3500 

17.3 

M7 III 

-0.3 

17 000 

15.9 

6300 

18.1 

F2 V 

3.6 

4500 

16.7 

2700 

17.8 







F8 V 

4.0 

4000 

17.0 

2500 

18.1 

B0 lb 

-6.1 

33 000 

11.5 

9500 

15.4 

G2 V 

4.7 

3500 

17.2 

2200 

18.2 







K3 V 

6.65 

2400 

18.4 

1700 

19.1 

WD 

8.0 

1500 

18.9 

1200 

19.2 

MOV 

8.8 

1500 

19.7 

1200 

20.0 







M8 V 

13.5 

500 

21.8 

450 

22.1 








errors) for a sample of main sequence stars selected from a small 
volume element in the Galaxy. Thus, main sequence fitting to such 
stellar samples, selected by means of the Gaia parallaxes, could 
serve as the basis for a 3D extinction map jKniidefooj) . The chem¬ 
ical composition of the stars may be used to refine the definition of 
an appropriate main sequence, and log g to refine the selection of an 
appropriate stellar sample. This approach has the further advantage 
that the C1B passbands may be used in addition to the C1M ones, 
and that stars along the entire main sequence will contribute to the 
estimation of extinction, including those located on the lower main 
sequence where extinction determination is otherwise difficult. 

The combination of Gaia and external IR photometric data 
provides independent determinations of extinction. During the lu¬ 
minosity calibration of 2MASS photometry for A9-G5 main se¬ 
quence stars. lKnude & Fabricius ( 2003 ) noticed significant redden¬ 
ing vectors in the Mj vs. (J — H) a b s diagram where Mj was 
estimated from Hipparcos parallaxes ty > 8.0 milliarcsec and 
a n /-7T < 0.11, assuming that reddening was negligible. These vec¬ 
tors are most pronounced for the A type stars but are also present 
for cooler ones. From these vectors the extinction law may be esti¬ 
mated even for various spectral classes if the stellar density is large 
enough. Similar diagrams may be constructed from Gaia astrom¬ 
etry and BBP and MBP for the optical region and for the near- 
infrared from the 2MASS, UKIDSS and VISTA surveys so that 
regions with abnormal extinction can be identified. 


8.2 Absolute luminosity and gravity 

Photometry will be crucial for absolute luminosity (or gravity) de¬ 
termination when relative parallax errors are larger than 10-20%. 
The distances at which the relative parallax error is ~10% are listed 
in Table[5]as a function of spectral type and luminosity class in the 
absence of interstellar extinction, i.e. similar to the Galactic pole di¬ 
rection, and in the case of an average Galactic plane interstellar ex¬ 
tinction of 0.7 mag kpc -1 . T he parallax a ccuracies are from Table 
8.4 in the Gaia Study Report fesAl200d) . Based on the Besan§on 
Galaxy model iRobin et all2003i) ~some 100-200 million stars are 
predicted to have parallaxes measured at the 10% accuracy level. 

Absolute magnitude calibrations can be established using the 
stars with accurate parallaxes and interstellar extinction, and ap¬ 
plied to more distant stars in the classical way, assuming that these 
are intrinsically similar to the accurate-parallax stars, which may 
not always be the case. Thus luminosities with precision om v ~ 


0.2-0.4 mag (or ai osg ~ 0.1-0.2 dex) from photometry are desir¬ 
able to match the uncertainties of the well measured parallaxes. 

The main contributors to the photometric luminosity deter¬ 
mination are the measurements of Balrner and Paschen jumps, 
the Fla index, the C1M395-C1M410 vs. W(CaT*) diagram and 
C1M515 combined with contiguous pseudo-continuum passbands, 
as already noted in Section[5] 

Figures fl2l 1131 and 1151 show that ai ogg ~ 0.1-0.2 dex is 
achievable for giants and early type stars at large distances and 
therefore that the study of Galactic structure on large scales (warp, 
spiral arms, outer halo, etc) is feasible. For nearby stars the pho¬ 
tometric determination of log g is also possible but of less interest 
because of the good parallax precision. For distant late type dwarf 
stars, i.e. with parallax errors larger than 10%, the precision of the 
photometrically determined log g will be not as good as for giants 
and early type stars, as expected. 

8.3 Chemical abundances 

High resolution, high SNR spectroscopic observations allow 
the most reliable chemical abundances determinations (e.g., 
fCavrel de Strobel~et lifcooib . RVS spectra will be used to deter¬ 
mine the atmospheric parameters, when possible. [Katz et alJ. t) 2004 ) 
estimate <T[m/h] t0 be smaller than 0.2-0.3 dex for stars with 
V ~ 14 and to improve by combining the spectra with photometric 
and astrometric data. Thus, atmospheric parameters for 10-25 mil¬ 
lion stars to V ~ 14-15 will be determined including individual 
abundances for 2-5 million stars to V ~ 12-13 (mainly Fe, Ca, 
Mg and Si for F-G-K stars ; N in hotter stars, such A-type stars; 
and information on C, N or TiO abundances for cool K and M-type 
stars). 

Although 10-25 million stars with RVS chemical abundance 
determinations may seem little compared to the 1 billion objects 
observed by Gaia , the spectroscopic chemical compositions are 
crucial as they will serve to calibrate the photometric data. This 
is important because for the vast majority of stars, chemical abun¬ 
dances will be derived from photometric data exclusively. 

The chemical abundances, and especially the relation between 
[M/H] and [a/Fe], are indicative of the initial star formation rate 
(SFR) and provide the rate of chemical enrichment of the Galaxy. 
Calculations by Maede r 1 1999 i show that the a-element abundance 
ratio starts to decrease a t high [M/H] a s the initial SFR increases. 
In agreement with this, INissen ' t l 1999 9 argues that the thick disk 
stars underwent chemical evolution with a high initial SFR such 
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that a relatively high metallicity ([M/H] ~ —0.4) was reached 
before the contribution by iron due to the SNe la decreased the 
[a/Fe] value. He also argues that the initial SFR in the thin disk 
has been slower and that the relative a-element to iron abundances 
started to decrease at lower metallicity ([M/H] ~ —0.6). Most 
of the halo stars were formed at an even lower initial SFR, so the 
decrease of [a/Fe] occurs at low metallicity ([M/H] ~ —1.2). 
Some of the halo stars have [a/Fe] abundances typical of the thick 
disk which points to a dual model for the halo formation: an inner 
part which had a high SFR and an outer part which experienced a 
slower evolution or was accreted from dwarf galaxies. 

Figures fl2lll3l and 1 151 show that Gaia photometry is able to 
match the expected spectroscopic precision up to about 1-2 kpc 
from the Sun depending on the Galactic direction (i.e. the redden¬ 
ing) for F-M dwarfs and subdwarfs. For giants in the red clump 
and the red giant branch, both being brighter, the same precision is 
attainable up to about 4, 7 and 12 kpc in the centre, anti-centre and 
orthogonal Galactic directions, respectively. This implies that the 
following issues (and other goals) can be addressed: the determina¬ 
tion of the Galactic chemical abundance gradient, the classification 
of the stars into different stellar populations, the characterization of 
halo streams, and the determination of the distance scale through 
the metallicity determination for RR-Lyrae. 

As expected, low metallicities are determined with less pre¬ 
cision than solar or higher metallicities. Extensive spectroscopic 
follow-up efforts from the ground will be necessary to determine 
the metallicities of stars with [M/H] < —3.0 with a sufficient 
accuracy to select the most interesting targets for high-resolution 
spectroscopy (Christlieb, private communication). 

Photometric chemical composition determinations are not free 
of difficulties. For instance, cool unresolved binaries tend to mimic 
a single star with a lower metal content. As an example, for G and 
K dwarfs with companions for which Am is 1.5-2 mag the sys¬ 
tematic metallicity error is A[M/H] ~ —0.4. This bias decreases 
to A[M/H] ~ —0.2 when Am = 3. 

As discussed in Sec. Q the first trials to determine [a/Fe] 
using the NextGen2 library ([M/H] € [—2, +0] and [a/Fe] € 
[—0.2, +0.8]) show that precision of 0.1-0.2 dex for G ~ 15-16.5 
can be achieved. The determination may be improved with the in¬ 
clusion of parallax information, because it constrains the value of 
log g, which may break the degeneracy in the variation of the Mg 
lb triplet due to luminosity and Mg abundance changes. For stars 
fainter than G ~ 16-17, deriving a-elements abundances from 
photometry is unlikely. 

8.4 Age 

The location of a star in the HR-diagram does not allow for a unique 
age determination as several combinations of chemical compo- 
sition ([M/H], [q/Fe], Y , etc) and age are possible (see e.g., 
iBinnev & Merrifieldil 1998^1 . Individually accurate ages require ex¬ 
tremely accurate determinations of luminosity, temperature, red¬ 
dening and chemical composition, which Gaia will provide. 

For ages around 10-14 Gyr, and assuming a given chemical 
composition, a change of log T e ff by ~ 0.01 dex for the stars in 
the turn-off translates to an age variation of 2 Gyr. Thus, tempera¬ 
tures must be known with precisions better than a few hundredths 
in log T e fj, i.e. 1.5-2% in T e ft, for the turn-off stars to determine in¬ 
dividually accurate ages. A variation of 0.3 dex in [M/H] is equiv¬ 
alent to a change of log T e g by ~ 0.01 dex for a given age. The 
variation of the a-elements abundance has a smaller impact on the 
age determination than the variation of [M/H]. A variation of +0.3 


dex in [a/Fe] leads to a variation of —0.006 dex in the temperature 
turn-off, yielding an uncertainty of about 1 Gyr in the age. 

In summary, an uncertainty of about 4 — 5 Gyr is estimated for 
the individual ages of the stars at the turn-off for the halo, thick and 
old-thin disk stars up to about 2, 3 and 5 kpc in the Galactic centre, 
anti-centre and orthogonal directions, respectively, as deduced from 
the PS performances in SectionQ 

The age determination for F-G sub-giant stars is quite insensi¬ 
tive to uncertainties in T e fj because the isochrones are almost hori¬ 
zontal. An uncertainty in Mv of about 0.15 translates to an uncer¬ 
tainty in the age of about 2 Gyr. Assuming additional uncertainties 
in [M/H] and [a/Fe] determinations of about 0.3 dex, the final 
prec ision of the individual ages is about 3 — 4 Gyr. 

iKucinskas et al. i ( 2003 ) showed that early-AGB stars can pro¬ 
vide ages as accurate as the turn-off stars for [M/H] > —1.5 when 
ciog T eff 0.01, criog g ^ 0.2, <T[m/h] ^ 0.2 and (Jeb_y ^ 0.03, 
thus allowing the age determination to larger distances than for 
turn-off and sub-giant stars. The Gaia capabilities in the case of 
metal-poor stars have not yet been investigated. 

Subsets of each Galactic population (such as globular clusters, 
open clusters, OB associations, a given halo stream, an identified 
merger, a moving group, etc.) can be treated statistically and mean 
ages and chemical abundances of the group can be obtained with 
much better precisions than those of the individual members. 

The prospects for the determination of the age-metallicity re¬ 
lation and the star formation history with Gaia are discussed by 
lHavwoo5 i2005ti . Assuming the errors on the astrophysical param¬ 
eters from Section^ the author recovers the simulated age metal¬ 
licity relation in the case of moderate extinction and for a slowly 
varying star formation rate. The details of the star formation his¬ 
tory, such as a series of prominent short bursts, are not recovered. 

8.5 QSOs 

QSOs play an important role in the Gaia mission as they will be 
used to construct the astrometric reference frame and they are of 
course interesting in their own right. However, since the QSO pop¬ 
ulation only represents 0.05% of the stellar population, building a 
secure QSO sample with no stellar contamination requires a very 
efficient rejection algorithm. Although proper motion, parallax and 
variability information will help in rejecting stars, these will be 
available with the required precision only at the end of the mis¬ 
sion. Therefore it is important to check the capability to classify 
QSOs using only photometric data. 

The QSO classification efficiencies of Supervised NNs and of 
MDMs (see Sec. l7.2> have been compared using synthetic data gen¬ 
erated from the BaSeL 2.2 library for stars ( ILeieune et alJI 19981) . 
from pure-hydrogen atmosphere models for white dwarfs (Koester, 
private communication) and from a library of QSO synthetic spec¬ 
tra Iciaeskens et afl l2005albl) . Properly trained NNs are found to 
be capable of rejecting virtually all stars, including white dwarfs. 
This is at the expense of the completeness level of the QSO sample, 
being only ~ 20% at G = 20. However, this provides a sufficient 
number of objects in the context of the non-rotating extragalactic 
reference system determination. MDMs provide a higher complete¬ 
ness level (~ 60% at G = 20), but with a correspondingly higher 
stellar contamination rate of the QSO sample. MDMs may be pre¬ 
ferred at high Galactic latitudes. Not surprisingly, reddened QSOs 
and weak emission-line objects are preferentially lost while high 
red-shift objects are most easily recognized. 

Assuming an object is a QSO, it is also possible to infer its 
red-shift from its photometric signature in the C1B+C1M pass- 
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bands. Unfortunately, there is a colour degeneracy in the QSO spec¬ 
tra which limits the expected precision of the technique to about 
| |Median — 0.2 in the range 0.5 < z spe c < 2 (computa¬ 

tions were done with a previous version of the Gaia PS but are 
not strongly sensitive to the adopted filter sets). 

More details on the identification and characterization of 
QSOs with Gaia can be found in lClaeskens et ah (, 2005 a). 


9 CONCLUSIONS 

To fully achieve the scientific goals of the Gaia mission it is es¬ 
sential to complement the astrometric and radial velocity measure¬ 
ments with highly accurate multi-colour photometry. This is nec¬ 
essary both for obtaining accurate astrometry and for a proper sci¬ 
entific interpretation of the stereoscopic census of the Galaxy Gaia 
will provide. The latter goal requires that the photometric system 
for Gaia is capable of astrophysically parametrizing 1 billion ob¬ 
jects across the entire Hertzsprung-Russell diagram in the presence 
of varying degrees of reddening and for the full range of chemi¬ 
cal compositions and ages of stars populating the Galaxy, and, of 
identifying peculiar objects and QSOs. 

There are no existing photometric systems capable of handling 
this very demanding task. We therefore set out to design a new sys¬ 
tem for the Gaia mission. The broad passbands implemented in the 
Astro instrument should fulfil the chromaticity requirements from 
the astrometric data processing while at the same time allowing for 
the astrophysical characterization of stars in very crowded regions. 
The main classification and astrophysical parametrization task is 
carried out with the medium passbands implemented in the Spectro 
instrument. 

The novel development we introduce in this work is the use of 
an objective ‘figure of merit’ (FoM) to compare different proposals 
for the Gaia photometric system. This FoM is based on the pre¬ 
dicted errors that can be obtained for the parametrization of stars 
in terms of T e g, log g, Av, [M/H], and [a/Fe]. These errors are 
calculated by using synthetic spectra to evaluate the sensitivity of 
each photometric passband to the astrophysical parameters. At the 
same time the degree of local degeneracy in the astrophysical pa¬ 
rameters is taken into account in the FoM. The overall FoM for a 
photometric system is calculated taking into account the priorities 
of the different scientific targets that were deemed to be most im¬ 
portant for addressing the core science case for Gaia. Using the 
FoM allows us to choose objectively which of the many proposed 
photometric systems is best in terms of reaching the astrophysical 
parameters error goals while at the same time having the least num¬ 
ber of local degeneracies. 

The calculation of the predicted posterior errors and the FoM 
as outlined in Appendix[A]can be applied to any photometric sys¬ 
tem and can also be used to predict parametrization errors that can 
be achieved with spectra of stars. In addition our method can be ex¬ 
tended to the optimization of systems that can be used for the study 
of high red-shift galaxies or quasars (for example to derive photo¬ 
metric red-shifts), as long as the dependence of the relevant SEDs 
on parameters such as red-shift or star formation rate are known. 

The main contribution of this work is the photometric system 
itself. This system has been designed based on our astrophysical 
knowledge and experience with ground-based systems and resulted 
in five (C1B) and fourteen (C1M) passbands implemented in the 
Astro and Spectro instruments, respectively. The number of pass- 
bands reflects the variety of targets (all types of stars, QSOs, galax¬ 
ies, solar-system objects), the required stellar astrophysical param¬ 


eters (T e ff, logp, Av, [M/H], [a/Fe]) and the need of breaking 
degeneracies (changes in two or more parameters may translate to 
the same changes in some spectral features but different changes 
in other ones). In summary, three of the broad C1B passbands are 
located to the left of H/3 line, on the Ha line and to the right of 
Paschen jump. The other two passbands fill the gaps so that full 
coverage of the whole spectral range of Gaia observations is pro¬ 
vided. A broad passband in C1M, implemented in the Spectro in¬ 
strument, provides the measurement of the UV flux at wavelengths 
blueward of the Balmer jump. Seven medium passbands are placed 
on the crowding of Balmer series in early-type stars, on the Call H 
line, on Mg I triplet and MgH band, on the Ha line, on one of TiO 
absorption bands in cool stars, and on the strong CN band in R- and 
N-type stars. An additional passband serves to measure the flux in 
the wavelength range covered by RVS. Finally, five passbands are 
devoted to the measurement of the pseudo-continuum. 

The unfiltered light measured in the Astro and Spectro in¬ 
struments provides photometric data in two very broad passbands 
(from 400 to 1000 nm and from 350 to 1025 nm). The G magnitude 
from Astro yields the highest SNR among all Gaia magnitudes, and 
hence is the most suitable for variability analysis. 

End-of-mission magnitude precisions have been estimated fol¬ 
lowing an ‘aperture photometry’ approach. A precision of 0.01 mag 
is obtained at V ~ 18 and 16 for the C1B and C1M passbands, 
respectively. In the case of the G passband and for stars brighter 
than V ~16 and ~14 measured with C1B and C1M, the pho¬ 
tometric precision is limited by how accurately the data are cal¬ 
ibrated. Although there is no detailed calibration model, a rough 
estimation of the number of involved calibration parameters and 
the number of available stars for calibration purposes shows that 
sub-millimagnitude photometric precisions are potentially achiev¬ 
able. A PSF-fitting approach shows that it is possible to deal with 
stellar densities up to about 2^1 x 10 5 stars/sq deg to magnitude 20 
with Spectro, if accurate positional information from astrometry is 
used. 

The performance of the Gaia photometric system with respect 
to the determination of astrophysical parameters for single stars has 
been evaluated using the ‘posterior’ errors from the FoM formalism 
and using algorithms designed for the astrophysical parametriza¬ 
tion of stars (based on MDM and NN). The differences and limita¬ 
tions have been discussed. Both approaches demonstrate that preci¬ 
sion goals in Section l4~3l are generally achieved to G ~ 17—18 mag 
and that the precision for each parameter depends on temperature, 
luminosity, chemical composition and interstellar extinction, and 
apparent magnitude, as expected. Reasonable uncertainties of 0.5- 
1 dex are obtained even for the chemical composition of metal-poor 
stars ([M/H]= —4.0). The determination of [a/Fe] to a precision 
of about 0.1-0.2 dex is possible down to G~ 16-17 for low to in¬ 
termediate temperature stars. Assuming a standard extinction law, 
the C1B and C1M passbands allow the determination of individual 
extinctions and hence the determination of T e fj is reliable, even in 
the presence of significant extinction, down to at least G = 18. By 
combining photometry and astrometry, a statistical approach to the 
construction of a 3D Galactic extinction map is also possible. Gaia 
will yield parallaxes with relative uncertainties lower than 10% for 
some 100-200 million stars allowing absolute magnitude determi¬ 
nations across the entire HR-diagram if the extinction is known. 
The C1B+C1M system is designed to provide the same luminos¬ 
ity errors (0.1-0.2 dex in log<?) for giants and early type stars at 
large distances ensuring the study of the Galactic structure on large 
scales. Individual ages of stars at the turn-off of the halo and old- 
thin and thick disk at distances up to about 2, 3 and 5 kpc in the 
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Galactic centre, anti-centre and orthogonal directions, respectively, 
can be determined with a precision of about 4-5 Gyr. The same 
precisions can be obtained for early-AGB stars, thus allowing the 
probing of larger distances than with turn-off and subgiant stars. 
For subgiants, our estimations yield slightly lower uncertainties of 
3-4 Gyr. Note that we quote uncertainties for individual ages. For 
a group of stars the mean age will be determined to much better 
precision. 

An additional merit of the Gaia photometric system is its abil¬ 
ity to discriminate QSOs from stars and white dwarfs. This will 
results in a sufficient number of objects for the definition of the 
non-rotating extragalactic reference system for proper motions. 

Finally, the performance of the C1B passbands when deal¬ 
ing with the chromaticity residuals has been evaluated. Our esti¬ 
mates show that the chromatic contribution to the parallax errors 
is 0.14-1.4 microarcsec. As these numbers are based on a worst- 
case assumption for the WFE and a somewhat simplistic calibra¬ 
tion model, a residual contribution of 1 microarcsec for stars and of 
few microarcsec for the QSOs is likely achievable. 

In summary, the photometric system C1B+C1M developed by 
the Photometry Working Group satisfies the mission requirements 
and was therefore adopted as the baseline by the Gaia Science 
Team. It has since been adopted by the ESA Project Team as the 
basis for the formal mission requirements to the industrial teams 
participating in ESA’s Invitation to Tender for Gaia. 

The experience from the development of the photometric sys¬ 
tem, including selecting the scientific targets and the ‘figure of 
merit' approach, puts us in a good position to rapidly optimize new 
payload proposals that will come from the selected industrial con¬ 
tractor in the course of 2006. 


ACKNOWLEDGMENTS 

We would like to acknowledge the contributions by B. Edvards- 
son, R. Lazauskaite, T. Lejeune and J. Sudzius and the collabora¬ 
tive effort of members of the Photometry Working Group within the 
Gaia mission preparation framework. We also thank P. Hauschildt, 
I. Brott and F. Castelli for making available to us SED libraries 
with a-elements enhanced abundances. The authors acknowledge 
the funding support from the following countries: Belgium, Den¬ 
mark (through the Danish Ministry of Science), Estonia (through 
the ESF grants No. 5003 and 6106), France, Germany (through 
the DLR and the Emmy Noether Programme of the DFG), Greece 
(through GSRT), Italy, Lithuania (through the Ministry of Edu¬ 
cation and Science), The Netherlands (through the Netherlands 
Research School for Astronomy NOVA and The Netherlands Or¬ 
ganisation for Scientific Research NWO), Russia, Spain (through 
MCYT under contract PNE2003-04352), Sweden, Switzerland and 
United Kingdom. For Denmark, Estonia, Lithuania and Sweden 
support from the Nordic Research Board [Ref. NB00-N030 and 
2001-2783-4] is gratefully acknowledged. A. Kucinskas acknowl¬ 
edges the Wenner-Gren Fellowship. Finally, we would like to thank 
the referee, M.S. Bessell, for a rapid response and constructive 
comments that helped improve the paper. 


REFERENCES 

The Gaia technical reports are available from the Gaia Photome¬ 
try Working Group website: 


http://gaia.am.ub.es/PWG/documents/MNRAS/ 

Arenou F., Babusiaux C., Chereau F., Mignot S., 2005, in Pro¬ 
ceedings of the Symposium ‘The Three Dimensional Universe 
with Gaia’, ESA SP-576, p. 335 

Bailer-Jones C. A. L., 1998, Technical report, Statnet - a feedfor¬ 
ward interpolation neural network, 
http://www.mpia.de/homes/calj/statnet.html 
Bailer-Jones C. A. L„ 2000, A&A 357, 197 
Bailer-Jones C. A. L., 2003, in Munari U., ed., ASP Conf. Ser. 
Vol. 298, GAIA Spectroscopy, Science and Technology, Astro¬ 
nomical Society of the Pacific, San Francisco, p. 199 
Bailer-Jones C. A. L„ 2004, A&A, 419, 385 
Bailer-Jones C. A. L., 2005, in Proceedings of the Symposium 
‘The Three Dimensional Universe with Gaia’, ESA SP-576, p. 
393 

Bessell M. S„ Brett J. M„ Scholtz M„ Wood P. R„ 1989, A&AS, 
77, 1 

Bessell M. S., Castelli F„ Plez B„ 1998, A&A, 333, 231 (erratum 
A&A, 337. 321) 

Bessell M. S„ 2005, ARA&A, 43, 293 

Binney J., Merrifield M., 1998, Galactic Astronomy, Princeton 
University Press 

Bridzius A., Vansevicius V., 2002, Ap&SS 280, 41 
Brott I., Hauschildt P.H., 2005, in Proceedings of the Symposium 
‘The Three Dimensional Universe with Gaia’, ESA SP-576, p. 
393 

Brown A. G. A., 2003, Results of the second cycle of blind testing, 
Gaia technical report ICAP-AB-004 
Brown A. G. A., 2005, in Proceedings of the Symposium ‘The 
Three Dimensional Universe with Gaia’, ESA SP-576, p. 377 
Brown A. G. A., Jordi C., Knude J., H0g E., 2004 Procedure for 
Photometric Systems Recommendation, Gaia technical report 
PWG-AB-003 

Cardelli J. A., Clayton G. C., Mathis J. S„ 1989, ApJ, 345, 245 
Cayrel de Strobel G., Soubiran C„ Ralite N„ 2001, A&A, 373, 
159 

Cellino A., Delbo M., DelTOro A., Zappala V., 2005, in Proceed¬ 
ings of the Symposium ‘The Three Dimensional Universe with 
Gaia’, ESA SP-576, p. 231 

Claeskens J.-F., Smette A., Surdej J., 2005a, in Proceedings of the 
Symposium ‘The Three Dimensional Universe with Gaia’, ESA 
SP-576, p. 667 

Claeskens J.-F., Smette A., Vandenbulcke L., Surdej J., 2005b, 
MNRAS, accepted 

Dollet C., Bijaoui A., Mignard F., 2005, in Proceedings of the 
Symposium ‘The Three Dimensional Universe with Gaia’, ESA 
SP-576, p. 429 

Drimmel R., Cabrera-Lavers A., Lopez-Corredoira M., 2003, 
A&A, 409, 205 

Drimmel R., Bucciarelli B., Lattanzi M. G., Spagna A., Jordi C., 
Robin A. C., Reyle C., Luri X., 2005, in Proceedings of the Sym¬ 
posium ‘The Three Dimensional Universe with Gaia', ESA SP- 
576, p. 163 

ESA-SCI(2000)4, 2000, GAIA: Composition, Formation and 
Evolution of the Galaxy 

Evans D. W., 2004, MBP Profile Fitting in Crowded Regions, 
Gaia technical report PWG-DWE-001 
Eyer L., 2005, in Proceedings of the Symposium ‘The Three Di¬ 
mensional Universe with Gaia’, ESA SP-576, p. 513 
Eyer L„ Mignard F„ 2005, MNRAS, 361, 1136 
Fiorucci M., Munari U., 2003, A&A, 401, 781 


24 Jordi et al 


FitzpatrickE. L„ 1999, PASP, 111, 63 

Fluks M. A., Plez B., The P. S., de Winter D., Westerlund B. E., 
Steenman H. C., 1994, A&AS, 105, 311 
Grenon M., Jordi C., Figueras F.. Torra J., 1999, An intermediate 
band system for GAIA, Gaia technical report MGUB-PWG- 
002 

Gunn F. E„ Stryker L. F„ 1983, ApJS, 52, 121 
Gustafsson B., Edvardsson B., Eriksson K., Jprgensen U.G., 
Mizuno-Wiedner M., Plez B., 2003, in Proceedings of the IAU 
Symposium 210, Piskunov N„ Weiss W.W., Gray D.F., eds., p. 
A4 

Haywood, M„ 2005, in Proceedings of the Symposium ‘The 
Three Dimensional Universe with Gaia’, ESA SP-576, p. 521 
Hastie T., Tibshirani R., Friedman J., 2001, The Elements of Sta¬ 
tistical Learning, Springer-Verlag, New York 
Hauschildt P. H„ Allard F„ Baron E„ 1999, ApJ, 512, 377 
Hauschildt P. H., Allard F., Baron E., Aufdenberg J., Schweitxer 
A., 2003, in Munari U„ ed„ ASP Conf. Ser. Vol. 298, GAIA 
Spectroscopy, Science and Technology, Astronomical Society of 
the Pacific, San Francisco, p. 179 
Heiter U., Piskunov N., Gustafsson B., Jordi C., Carrasco J. M., 
2005, in Proceedings ‘13th Cambridge Workshop on Cool Stars, 
Stellar Systems and the Sun’, ESA SP-560. p. 635 
H0g, E. 1993, Astrometry and Photometry of 400 Million Stars 
Brighter than 18 mag, in IAU Syrnp. 156. 37-45, Kluwer, 
Netherlands 

H0g, E. 1998, Spectrophotometry or Filter Photometry for Gaia?, 
Gaia technical report SAG-CUO-045 
H0g E., 2005, Summary of sampling schemes for ASM, AF, BBP, 
SSM and MBP, Gaia technical report GAIA-CUO-151V3 
H0g E., Knude J., 2004a, Two new MBP systems: K2M and K3M, 
Gaia technical report GAIA-CUO-152 
H0g E., Knude J., 2004b, Four new MBP systems: K4M, K5M, 
K6M, K7M, Gaia technical report GAIA-CUO-155 
Jordi C., Carrasco J. M., 2004a, BBP systems proposals, Gaia 
technical report UB-PWG-026 

Jordi C., Carrasco J. M„ 2004b, MBP systems proposals, Gaia 
technical report UB-PWG-027 

Jordi C., Grenon M., Figueras F., Torra J., Carrasco J. M., 2003, 
The 2F Geneva-Barcelona medium band photometric system for 
GAIA, Gaia technical report UB-PWG-011 
Jordi C., Carrasco J. M., H0g E., Brown A. G. A., Knude J., 
2004a, BBP photometric systems evaluation, Gaia technical re¬ 
port UB-PWG-028 

Jordi C., Carrasco J. M., H0g E., Brown A. G. A., Knude J., 
2004b, BBP photometric systems evaluation, Gaia technical re¬ 
port UB-PWG-029 

Jordi C., Figueras F., Torra J.. Carrasco J. M.. 2004c, The 3F 
Geneva-Barcelona medium band photometric system for Gaia, 
UB-PWG-016 

Jordi C., Figueras F., Carrasco J. M„ Knude J., 2004d, Gaia Sci¬ 
entific targets for PS design, Gaia technical report UB-PWG- 
009v 1.1 

Jordi C., Knude J., Carrasco J. M., Figueras F., 2004e, Gaia Sci¬ 
entific targets for PS design: quantification of priorities, Gaia 
technical report UB-PWG-015vl.2 
Kaernpf T.A., Willemsen P.G., Bailer-Jones C.A.L., 2005, in Pro¬ 
ceedings of the Symposium 'The Three Dimensional Universe 
with Gaia’, ESA SP-576, p. 441 
Kaltcheva N„ Knude J., Georgiev V., 2003, A&A, 407, 377 
Katz D., Soubiran C., Cayrel R., Adda M., Cautain R., 1998, 
A&A, 338, 151 


Katz D„ et al., 2004, MNRAS, 354, 1223 
Knude J., 2002, Ap&SS, 280, 97 
Knude J., Fabricius C., 2003, Baltic Astronomy, 12, 508 
Knude J., H0g E., 2004, A new MBP system: KIM, Gaia techni¬ 
cal report GAIA-CUO-148 

Kolka I., Eenmae T., Hirv A., Tuvikene T., Kama M., 2005, in 
Proceedings of the Symposium ‘The Three Dimensional Uni¬ 
verse with Gaia', ESA SP-576, p. 543 
Kucinskas A., Lindegren L., Tanabe T., Vansevicius V., 2003, in 
Munari U., ed., ASP Conf. Ser. Vol. 298. GAIA Spectroscopy, 
Science and Technology, Astronomical Society of the Pacific, 
San Francisco, p. 415 

Kucinskas A., Hauschildt P. H.. Ludwig H.-G., Brott I., Van¬ 
sevicius V., Lindegren L., Tanabe T., Allard F., 2005, A&A 442, 
281 

Kurucz R. L„ 1979, ApJS, 40, 1 
Lejeune T., Cuisinier F., Buser R., 1998, A&AS, 130, 65 
Lindegren L., 2001, GAIA chromaticity calibration and design 
of the Broad Band Photometry system, Gaia technical report 
GAIA-LL-039 

Lindegren L., 2003a, Photometric systems for Gaia’s Broad band 
Photometer, Gaia technical report GAIA-LL-045 
Lindegren L., 2003b, Optimizing Gaia’s Photometric System: 
Thoughts on distance measure and figure of merit, Gaia tech¬ 
nical report GAIA-LL-047 

Lindegren L., 2003c, Gaia chromaticity calibration and the BBP 
filter shape, Gaia technical report GAIA-LL-049 
Lindegren L., 2004, Use of parallax information in the Photomet¬ 
ric System design, Gaia technical report GAIA-LL-054 
Maeder A., 1999, Proceedings of ‘The Galactic halo: from glob¬ 
ular clusters to field stars’, A. Noels et al (eds.), Universite de 
Liege, p. 1 

Malyuto V., 2005, in Proceedings of the Symposium ‘The Three 
Dimensional Universe with Gaia’, ESA SP-576, p. 453 
Megessier C., 1995, A&A, 296, 771 

Mignard F.. 2005, in Proceedings of the Symposium ‘The Three 
Dimensional Universe with Gaia’, ESA SP-576, p. 5 
More D„ Munari U„ 2000, A&AS, 147, 361 
Munari U., 1999, Baltic Astronomy, 8, 123 
Munari U., Sordo R., Castelli F., Zwitter T., 2005, A&A, 442, 
1127 

Nissen P.E., 1999, Proceedings of ‘The Galactic halo: from glob¬ 
ular clusters to field stars’, A. Noels et al (eds.), Universite de 
Liege, p. 125 

Nurmi P„ 2005, Ap&SS, 298, 559 

Perryman M. A. C., de Boer K. S., Gilmore G., H0g E., Lattanzi 
M. G., Lindegren L., Luri X., Mignard F., Pace O., de Zeeuw 
P. T„ 2001, A&A, 369, 339 
Pickles A. J., 1998, PASP, 110, 863 

Press W. H., Teukolsky S. A., Vetterling W. T., Flannery B. P, 
1992, Numerical Recipes in Fortran: The Art of Scientific Com¬ 
puting, Second Edition, Cambridge University Press 
Robin A. C., Reyle C., Derriere S., Picaud, S., 2003, A&A 409, 
523 

Robin A. C., Reyle C„ Picaud S„ Schultheis M„ 2005, A&A 430, 
129 

Straizys, V. 1992, ‘Multicolor Stellar Photometry’, Pachart Pub¬ 
lishing House, Tucson, Arizona 

Straizys V., 2004, A revised version of the BBP system for Gaia: 

V1B, Gaia technical report GAIA-VILN-002 
Straizys V., H0g E., 1995, in Proceedings of the Joint RGO-ESA 
workshop ‘Future Possibilities for Astrometry in Space', ESA 



SP-379, 191 

Straizys V., H0g E., Vansevicius V., 2000, The 3G photometric 
system, Gaia technical report SAG-CUO-78 
Straizys V., Zdanavicius K., Lazauskaite R., 2004, A new ver¬ 
sion of the MBP system for Gaia: V2M, Gaia technical report 
G Al A-VILN-001 

Tautvaisiene G., Edvardsson B., 2002, Ap&SS, 280, 143 
Tautvaisiene G., Edvardsson B., 2005, in Proceedings of the Sym¬ 
posium ‘The Three Dimensional Universe with Gaia’, ESA SP- 
576, p. 549 

Tautvaisiene G., Edvardsson B., Bartasiute S., 2003, Baltic As¬ 
tronomy 12, 532 

Torra J., Chen B., Figueras F., Jordi C.. Luri X., 1999, Baltic As¬ 
tronomy, 8, 171 

Vansevicius V., Bridzius A., 2002, Gaia photometric perfor¬ 
mances, Gaia technical report GAIA-VIL-008 
Vansevicius V.. Bridzius A., 2003, in Munari U., ed., ASP Conf. 
Ser. Vol. 298, GAIA Spectroscopy, Science and Technology, As¬ 
tronomical Society of the Pacific, San Francisco, p. 41 
Wahlgren G. M., Lundqvist M., Kucinskas A., 2005, in Proceed¬ 
ings of the Symposium ‘The Three Dimensional Universe with 
Gaia’, ESA SP-576, p. 619 
Wilkinson M. I., et al„ 2005, MNRAS 359, 1306 
Willemsen P. G., Bailer-Jones C. A. L., 2005, Determination of a- 
abundances from Gaia photometric data, Gaia technical report 
ICAP-PW-008 

Willemsen P. G., Kaempf T. A., Bailer-Jones C. A. L., 2005a, 
Neural network parametrization performances for the C1M and 
C1B photometric systems, Gaia technical report ICAP-PW-006 
Willemsen P. G., Kaempf T. A., de Boer K. S., Bailer-Jones 
C. A. L., 2005b, Using distance information in the process of 
automated parametrization - preliminary results, Gaia technical 
report ICAP-PW-005 


APPENDIX A: FIGURE OF MERIT 

The definition and the procedure for calculating a figure of merit 
for a photometric system was proposed bv lLindegreri f' 2003bh and 
we provide the details here. 

The achievable errors <Tik, P ost for a particular ST i can be cal¬ 
culated using the so-called sensitivity matrix Si. The elements of 
this matrix are the partial derivatives d(f>ij/dpk, where (j>ij is the 
(noise-free) normalized flux in filter j (in photon counts) for ST i , 
and pk stands for the astrophysical parameter k. These derivatives 
describe how the flux in each filter changes in response to a change 
in astrophysical parameter (AP) k. Consider the AP determination 
as a linearized least squares estimation of A p, the improvement of 
the AP vector. The observation equation for ST i resulting from the 
flux measured in filter j reads: 

,'f lj Api + • • ■ + A pk = A (f>ij ± tij , (Al) 

dp i dpK 

where A <f>ij = 4>ij,ohs — (pij ( p ) is the difference between the 
observed and predicted flux and tij indicates the flux uncertainty. 
Observation equations of unit weight are formed through division 
by tij, whereupon normal equations are formed in the usual man¬ 
ner. The linearization is assumed to be made around the true param¬ 
eter vector p, so that the resulting update A p has zero expectation. 
Then, given the variance-covariance matrix Cj, = diag(e;.,) of the 
observed fluxes, the variance-covariance matrix of the estimated 
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AP-vector p post is given by the inverse of the normal equations 
matrix: 

C p , pos t = (S T C^S)~ 1 , (A2) 

where the matrices C^post, C$, and S are defined for each ST i 
separately. The diagonal elements [C' p , p0 st]fcfc = cr;fc !post of this 
matrix are the sought after achievable errors for a PS. In reality de¬ 
generacy among the APs will often make the matrix (S T Sf 1 
singular or near-singular, resulting in infinite or very large aik,post 
as computed from equation IA2> . This can be avoided by adding a 
suitable positive definite matrix B which makes the whole right- 
hand side positive definite; thus: 

C p , p OS t = (B + S T C ^ ) 1 Sf 1 . (A3) 

B is the a priori information matrix of the APs. In the absence of 
any other information on the APs we have B = diag(<T^. 2 p rior ). 
This matrix plays an important role and can be used to incorpo¬ 
rate prior information on the astrophysical parameters. This matrix 
can also be used to incorporate constraining information from the 
parallax measurements. When the photometric data does not pro¬ 
vide any relevant information on a given AP pk (either because the 
flux variances in C $ are too large or because the elements of the 
sensitivity matrix S are too small), then Oik, po st — Oik, prior- 

In practice the derivatives d(j>ij/dpk are calculated numeri¬ 
cally from simulated photometric data. The calculation thereof re¬ 
quires (synthetic) spectral energy distributions (SEDs) of the STs 
and a noise model for the photometric instruments. 

The figure of merit is now calculated as follows. For each ST 
i and astrophysical parameter k (Ay, [M/H], logg, T e fj, [a/Fe], 
...) the performance of the photometric system is measured by the 
ratio o ik, P ost /cAfcjgoai• The figure of merit Qi for each ST i is then 
defined as: 

Qi — 5 'ITk f {oik, post / Oik, goa\) , (A4) 

k 

where Wk indicates the relative weight of each astrophysical pa¬ 
rameter (with 52 & Wk = 1 ) and f(x) is a non-linear function of 
x = Oik ,post /<Xik ,goal with a break around 1. The function used is: 


/Or) = (1 + x 2n )~ 1/n . (A5) 


The global figure of merit (summed over all STs) is the weighted 
and normalized sum: 


52 ; 


t Vi 


(A6) 


where the weights Wi indicate the priority of each ST i. The value 
of Q indicates how close the performance of a photometric system 
is to being ‘ideal’ (when Q = 1). This global figure of merit is 
then calculated for each photometri c syste m using the error goals, 
weights and priorities defined in lJordi et al J i2004a ). 

The value of n in the function f(x) from equation 1\5> de¬ 
termines how much weight good performance (x ^ 1) gets in the 
FoM as opposed to bad performance (x > 1). We are looking for a 
PS that achieves the error goals over all of AP space and we want 
to avoid giving a high rank to a system that is extremely good in 
only one corner of AP space but very bad everywhere else. The lat¬ 
ter will be highly ranked for n = 1 but not for large n. However, 
for very large n (for which f(x) approaches a step function with 
f(x) = 0 for x > 1) we will not be able to tell the difference be¬ 
tween two PSs that are equally bad wherever x > 1 but with one 
PS being better than the other when x < 1. We chose the value of 
n = 3 as a good compromise between these two extremes. 
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We end with a few remarks about the sensitivity matrix. The 
columns of Si are the gradient vectors that describe the changes of 
the fluxes with respect to a change in an AP. Thus Si contains all 
the information needed to characterize the behaviour of the pho¬ 
tometric system in the data space (or filter flux space) near ST i. 
In the ideal case the gradient vectors would be aligned with the 
coordinate axes of the data space. That is each flux measurement 
<j>ij is sensitive to one (and only one) AP. In practice the gradient 
vectors are of course not aligned with the coordinate axes which 
means that each AP influences some linear combination of fluxes. 
This also means that even if all the gradient vectors are orthogo¬ 
nal the errors in the APs will still be correlated because any error 
on a measured flux <f>ij will influence multiple AP determinations. 
These correlations are correctly taken into account in the calcula¬ 
tion of the variance-covariance matrix of the estimated AP-vector 
and are reflected in non-zero off-diagonal elements in this matrix. 

A further complication that will occur in practice is that the 
gradient vectors will not be orthogonal to each other. This means 
that there will be degeneracies between the APs when we try to es¬ 
timate them. A well-know example is the degeneracy between T e fr 
and Av if only the continuum of the spectrum is measured. The 
behaviour of the PS depends on the noise and we should in fact 
consider the noise-weighted gradient vectors which have compo¬ 
nents 1 /tij x d(j)ij/dpk- When the orthogonality is defined with 
respect to the noise-weighted gradient vectors, non-orthogonal gra¬ 
dient vectors will lead to larger correlations between the errors in 
the estimated APs and they will also cause the standard errors on 
each AP to increase. Without going into the mathematical details 
this can be appreciated if one considers that any degeneracy be¬ 
tween two APs will make it more difficult to attribute flux changes 
to either of them thus increasing the uncertainty in both parameters. 

This means that a PS that contains larger degeneracies will 
get a lower figure of merit because of the increased aik, P ost for 
degenerate APs. Note, however, that the FoM used in this sense is 
a measure of how good a particular PS is at locally separating stars 
with different APs along orthogonal directions. The FoM does not 
take global degeneracies into account, where very different parts of 
AP space are mapped onto each other in filter flux space. 



